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To investigate heavy fermion behavior in the vanadium spinel LiV2 04, we start from a three- 
orbital Hubbard model on the pyrochlore lattice and derive its low-energy effective Hamiltonian by 
an approach of real-space renormalization group type. We first derive the effective Hamiltonian nu- 
merically, and then succeed in representing the results into an analytic form with physical operators 
for low-energy degrees of freedom in tetrahedron unit. The effective Hamiltonian is defined on the 
coarse grained lattice, i.e., face-centered cubic (f.c.c.) lattice , and it operates in a restricted Hilbert 
space defined in terms of a specific molecular orbital Ti in the unit. One important tetrahedron 
configuration has a three- fold orbital degeneracy and spin S = 1, and correspondingly, the effective 
Hamiltonian has spin and orbital exchange interactions of Kugel-Khomskii type as well as correlated 
electron hoppings. The coupling constants in the effective Hamiltonian are determined from the nu- 
merically obtained renormalized Hamiltonian, and also by means of perturbation. We calculate and 
analyze low-energy states of the effective Hamiltonian for the unit of four coupled tetrahedra both 
analytically and numerically. Effective hopping elements in the effective Hamiltonian are renormal- 
ized to about 1/10 of the original hopping integral. It is important that different virtual processes 
make opposite contributions to the exchange term, and consequently the coupling constant is given 
by a remaining small value. This is particularly prominent in the spin-spin channel, where ferromag- 
netic double exchange processes compete with antiferromagnetic superexchange processes. Another 
important point is that various spin and orbital exchange processes are competing to each other. 
Together with geometrical frustration of the effective f.c.c. lattice, these two features result in nearly 
degenerate three lowest-energy states of different types in the four coupled tetrahedra, and each of 
the three has a finite degeneracy in spin and/or orbital. We also calculate spatial correlations of spin 
and orbital and found that short-range spin-spin correlations are strongly entangled with orbital 
configurations. This indicates that a large remaining entropy at low temperature is related to slow 
coupled fluctuations of spin and orbital. These results suggest the absence of phase transition in 
spin and orbital spaces down to very low temperatures and their large fluctuations in the low-energy 
sector, which are key issues for understanding the heavy fermion behavior in L1V2O4. 

PACS numbers: 71.27.+a, 71.10.Fd 



I. INTRODUCTION 



The vanadium spinel LiV2 04 is the first heavy fermion 
compound discovered in d-electron systems 1 -. For about 
a decade, various experimental and theoretical efforts 
have been made to understand its heavy fermion be- 
haviors. Low-temperature properties such as specific 
heat, magnetic susceptibility, electrical resistivity and 
Hall coefficient seem to be explained by the quasiparticle 
picture with a large effective massi 2 -^! 5 - Corresponding 
to these low-temperature behaviors, the electronic spec- 
tral function develops a peak above the Fermi energy 
at low temperature observed in the laser photoemission 
spectroscopy^. All these low-temperature behaviors are 
characterized by one energy scale T* ~ 30 K. In con- 
trast to these low-temperature properties, L1V2O4 ex- 
hibits bad metallic behaviors at higher temperatures-. 
The temperature (T) dependence of magnetic properties 
is also interesting. The size of magnetic moment changes 
from a mixed- valent value (V 3+ +V 4+ ) to a smaller value 
at about 500 K£&& Neutron experiment showed that 
spin fluctuations x(q) change their spatial correlations 
at around T = T* . In the higher-temperature region, 



spin fluctuations show a broad peak at q = 0, and 
this is considered as a consequence of double exchange 
interactions. In the lower-temperature region, neutron 
experiments exhibited that the peak position shifts to 
|q| = |Q*| ~ 0.6 A^ 1 . 10 ' 11 Therefore, it is important 
to investigate the competition of ferro- and antiferro- 
magnetic fluctuations to understand the low-temperature 
heavy fermion behaviors. 

Regarding theories, it has been discussed that the 
heavy fermion behaviors originate from the Kondo 
effect ) 12 ' 13 : 14 the inter-orbital Coulomb interaction^ 
the spin-orbital fluctuations ) 16 ! 17 the frustrations ; 18 ! 19 
and dimensional crossover from coupled one-dimensional 
chains to three dimensions^ 

Band structure calculations show that the Fermi 
surfaces of L1V2O4 are composed by d-electrons t2 g 
orbitalsi 12 ' 21 ! 22 ' 23 ! 24 These t 2g orbitals split into ai g and 
e g orbitals due to the trigonal distortion of surrounding 
oxygen atoms. Anisimov et al. proposed the Kondo 
effect scenario^ 2 - that the <ii g orbital plays a role of lo- 
calized electron and interacts with conduction electrons 
in e g orbitals. A recent study of cluster dynamical mean 
field theory 25 claims that the heavy fermion behaviors are 
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related to the criticality of orbital-selective Mott transi- 
tion of dig electrons based on the analysis using a sim- 
plified two-orbital Hamiltonian. Yushankhai et al. ana- 
lyzed the low-temperature Q* spin fluctuations observed 
in the neutron scattering experiment by employing a phe- 
nomenological self-consistent renormalization theory of 
spin fluctuation 26 ' 27 . They succeeded in fitting the neu- 
tron data qualitatively, but understanding of the micro- 
scopic aspects of magnetic fluctuations and heavy quasi- 
particles are desired. Despite of these efforts, the compe- 
tition of ferro- and antiferro-magnetic interactions and 
crossover behaviors in the temperature dependence of 
susceptibility are not fully understood and it is desired to 
clarify how to describe the quasiparticles on a frustrated 
pyrochlore lattice and whether the frustration plays an 
important role for the realization of heavy fermion be- 
haviors. 

In this paper, we focus on the coupling of orbital de- 



calculation o 12 i 21 ' 22 i 23 ' 24 point out that the electronic den- 
sity of states near the Fermi energy consists mainly of the 
d-electron ti g orbitals on vanadium sites. In the spinel 
structure, the vanadium sites form a three-dimensional 
pyrochlore lattice and the unit cell contains four vana- 
dium atoms which form a tetrahedron as shown in Fig. 
[T] The electron hopping processes can be described by 
the effective V-V hoppings. Effects of V-0 hoppings are 
included as a renormalization of V-V hoppings. There is 
trigonal distortion in the lattice due to O ion displace- 
ment. This lifts three-fold degenerate i 2g orbitals into 
a\ g (singlet) and e g (doublet). The vanadium valence is 
y3.5+ j n avera g e a nd this corresponds to 1.5 electrons 
per atom, i.e., quarter filling of ti g orbital. 

The Hamiltonian we will investigate in this paper is a 
three-orbital £ 2g Hubbard model on the pyrochlore lattice 
with trigonal splittings, 
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starting with a microscopic model on the pyrochlore lat 
tice. We discuss its interplay with spin and charge de- 
grees of freedom and its spatial correlations beyond a 
tetrahedron cluster. To examine the role of orbital de- 
grees of freedom explicitly, we will use a three-orbital 
Hubbard model without assuming that a\ g orbital is lo- 
calized. Since the unit cell contains four vanadium atoms 
and each vanadium atom has three i 2g orbitals, straight- 
forward calculations are not applicable. In this paper, 
we shall employ a real space renormalization group ap- 
proach to extract a low-energy effective Hamiltonian for 
tetrahedron units^ 2 - The effective Hamiltonian is t- J-like 
model: localized spin-one and orbital-triplet degrees of 
freedom are coupled via exchange interactions and mo- 
bile electrons with three-fold orbital degeneracy hop be- 
tween tetrahedron units. Using this effective model, we 
will discuss low-energy electron itineracy and compet- 
ing interactions of spin and orbital degrees of freedom in 
LiV 2 4 . 

This paper is organized as follows. In Sec. HH we show 
the starting microscopic model used in this paper. Then 
in Sec. IIII1 we will demonstrate the results of the exact 
diagonalization to find low-energy degrees of freedom in 
one tetrahedron unit. In Sec. IIV1 we will discuss a pos- 
sible effective Hamiltonian which can describe the low- 
energy sector. In Sec. [V] we will show the four-unit di- 
agonalization results calculated by using the low-energy 
states in the one-tetrahedron calculations. In Sec. IVI1 
the low-energy physics is analyzed by the perturbative 
approach from the strong coupling limit. Finally, we dis- 
cuss the effective model relevant to L1V2O4 in Sec. IVIII 
and summarize the present paper in Sec. IVIIIl 



II. MODEL 

We start with describing a realistic microscopic 
model of electronic structure for the vanadium spinel 
LiV 2 04. In LiV 2 04, the first principle band 
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where d iaa is a d-electron creation operator with the or- 
bital a(= xy, yz or zx) and the spin cr(=t or |) at the site 
i, and its number operator is defined as n- lacl = d iaa di aa . 
The electron hoppings ty are limited to the nearest 
neighbor sites and fj, is the chemical potential. The 
trigonal splittings are included in t?f . For the inter- 
action term, we use standard onsite Coulomb interac- 
tions without pair hopping terms as in other studie o 16 ' 17 . 
The rotational symmetry of Coulomb interaction is sat- 
isfied by the relation U = U' + J and we will use U, 
U' and J satisfying this condition throughout this pa- 
per. We choose the nearest neighbor tight-binding pa- 
rameters ty by setting Slater-Koster parameters 3 ^ as 
tfj = (dda) = -0.527 eV, U = (ddn) = -0.085 ~ 
-0.13 eV, and t s = {ddS) = 0.25 eV for a-, tt- and 
(5-bond, respectively. We also fix the trigonal splitting 
A = e £g — e ais — 0.02 eV. Although A was estimated 
to be ~ 0.1 eV by the band calculation 12 , it does not 
directly correspond to the "microscopic" A we use^i 

Before examining the effects of electron correlations, 
we study the electronic structure of the non-interacting 
case. Each unit cell contains four vanadium atoms and 
each atom has three orbitals and therefore there are 
twelve bands in total. Their energy dispersion and the 
non-interacting density of states are shown in Fig. [2] 
for = —0.085 eV. Apart from the high-energy re- 
gions the overall features are in good agreement with the 
first-principle band structure calculation a 12 ' 21 ! 22 ' 23 . It is 
noted that the weight of a\ g orbitals is larger than the 
one of e g near the Fermi energy. Among six electrons 
per unit cell, the occupation numbers are n ai = 1.18 
and n e = 4.82 per tetrahedron in this parameter set. 

For later purpose, let us first consider molecular or- 
bitals in a single tetrahedron unit cell. Qualitatively, the 
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FIG. 1: (Color online) (a) The lattice structure of pyrochlore 
lattice. The unit cell contains four sites forming a tetrahe- 
dron, (b) Tetrahedron units in the pyrochlore lattice form an 
face-centered cubic (f.c.c.) lattice. 

energy levels of the molecular orbitals correspond to the 
band energies at the T point. There are twelve molecu- 
lar orbitals in total. The unit cell has the point group 
symmetry The twelve orbitals constitute five multi- 
plets labeled by irreducible representations of group, 
A 1: E, Ti, and T 2 (+) . Since there are two T 2 repre- 
sentations, we distinguish them by (— ) and (+). These 
orbitals are listed in Table U and the wavefunctions of 
the molecular orbitals are shown in Appendix [XJ We la- 
bel the irreducible representations by T and define the 
energy level as er- In this paper, we choose t n value 
such that the -level is higher than T^ - ^ -level and we 
show the t„ dependence of the energy at the T point in 
Fig. [31 We will see in Sec. IVII that the position of Al- 
and ^-levels are important for low-energy properties 
and this is sensitive to t n . The parameter region we dis- 
cuss in this paper is i OT = —0.085 ~ —0.13 eV and there 
£e < £rp(-) < ea x < £Ti < £r r (+). Large contribution of 
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ax g orbital near the Fermi surfaces comes from T% ' and 
A\ orbitals. 



III. MANY-ELECTRON EIGENSTATES OF 
ONE TETRAHEDRON 

In this section, we include Coulomb interactions in the 
Hamiltonian {1} and investigate its eigenstates in a single 
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FIG. 2: (Color online) Tight-binding dispersion and density 
of states calculated by the non-interacting Hamiltonian. The 
Fermi energy is indicated by dotted line. 
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FIG. 3: Non-interacting one-particle energies at Y point vs 

t 7T • 



tetrahedron. The results obtained in this section provide 
an insight to understand the high-temperature proper- 
ties of L1V2O4. Moreover, the many-body wavefunctions 
obtained in this section become good bases for the low- 
energy effective model which will be discussed in Sec. IIVI 
and four-tetrahedron calculations in Sec. [V] 

In the first part of this section, we will show the energy 
spectra of one tetrahedron unit calculated by exact diag- 
onalization for typical sets of parameters in the Hamilto- 
nian |T]). For L1V2O4, the average d-electron number is 
1.5 per vanadium site. This corresponds to six electrons 
per tetrahedron. We will discuss the energy spectra for 
the total d-electron number rid = 4, 5, 6 and 7. Then, in 



4 



TABLE I: Molecular orbitals of tetrahedron unit and their 
energy. D = [(§[i CT + U] + t s + f) 2 + 2(U - tg - f ) 2 ]i 
The third and the fourth columns represent the weight of ai g 
and e g atomic orbitals, respectively for = —0.085 eV. The 
values for two T2 representations depend on the hopping and 
trigonal splitting parameters but their sums are constant. 
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the second part, we will show the ground state phase dia- 
gram in the parameter space of Coulomb interactions. In 
the third part, we will also investigate how to construct 
low-energy spectra in terms of molecular orbital bases. 
In the final part of this section, we evaluate the temper- 
ature dependence of the thermodynamic quantities such 
as spin susceptibility, charge susceptibility and entropy. 



A. Energy spectra of one tetrahedron 

In this subsection, we show energy spectra of one tetra- 
hedron unit. The low-energy eigenstates of one tetrahe- 
dron unit will be used as bases of the discussion in later 
sections. We will also discuss thermodynamic properties 
such as spin susceptibility and entropy, and in order to 
examine the thermodynamic quantities in the whole tem- 
perature region, we need all the important eigenstates. 
Since the Hilbert space is very large, we restrict our- 
selves in the subspace for total d-electron number < 7 
in the tetrahedron unit and d°, d 1 , d 2 , and d 3 configura- 
tions on each vanadium atom. The other configurations 
d n with n > 4 are not taken into account. We have 
checked the validity of this truncated calculation at least 
for the purpose of discussing the low-energy properties 
of this system by comparing low-energy eigenvalues cal- 
culated by this truncated calculation with those of the 
full Lanczos method. The configurations included in the 
truncated calculations contain physical processes such as 
super- and double-exchange interactions which are im- 
portant when considering the low-energy properties of 
this system. The numerical diagonalizations were car- 
ried out with the open boundary condition, utilizing the 
spin rotational symmetry which reduces the maximum 
matrix size down to ~ 46000 for rid — 1 and the total 
spin 5=1/2. 

Figures S] (a)-(c) show the energy eigenvalues in each 
subspace of total spin S for U = 1.5 eV, U' = 1.3 eV, 
J = 0.2 eV, and t n = —0.085 eV. The numbers shown 
denote the degeneracy of each eigenstate that arises from 
the point group symmetry. The ground state of rid — 6 
is total spin 5 = 1 and orbital triplet. Figures 2] (d)-(f) 
show the energy eigenvalues for the same parameters in 
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FIG. 4: Energy spectra in each subspace of rid and spin (25*) 
for U = 1.5 eV, U = -0.085 eV. (a) n d = 5, (b) n d = 6, 

(c) n d = 7, (d) n d = 5, (e) n d = 6, and (f) n d = 7. The 
numbers denote the orbital degeneracy, (a)-(c): J = 0.2 eV, 

(d) -(e):J = 0.6 eV. 



(a)-(c) except J = 0.6 eV. As we increase J, the energies 
of the large spin states become lower. The ground state 
in rid — 6 sector is fully polarized. In the rid — 5 and 7 
sectors, however, the ground states are not fully polarized 
states. It is noted that the energy differences between 
different rid for J = 0.6 eV are smaller than those for 
J = 0.2 eV. This implies that charge fluctuations are 
enhanced as J is increased, which will be discussed in 
Sec. QTIDl 



B. Ground state phase diagram for one 
tetrahedron 



In Fig. we show the ground state phase diagram of 
the rid = 6 space in the U — J plane for — —0.085 
and —0.12 eV. Note that the region J > U is unphys- 
ical. These results are obtained for the full Hamilto- 
nian |T]) without truncating the Hilbert space. We use 
a usual Lanczos method to calculate the eigenenergies 
of the ground and the first excited states. There are five 
phases dependent on U, J and t w . Their total spin S and 
point group irreducible representation T are determined. 
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FIG. 5: Ground state phase diagram for one tetrahedron with 
n d = 6. U' = U-J. (a) U = -0.085 eV, and (b) -0.12 eV. 
The dotted lines are the boundary between different ground 
states. The large-J region in (a) (J > 1.0 eV) is similar to 
the corresponding region in (b) 
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FIG. 6: Schematic interpretation of ground states in the 
molecular orbital picture. 
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FIG. 7: Ground state phase diagram for rid = 5 and 7. V = 
U — J. (a) rid = 5 and = —0.085 eV. (b) rid = 5 and 
-0.12 eV. (c) n d = 7 and U = -0.085 eV. (d) n d = 7 and 
£rr = —0.12 eV. The dotted lines are the boundary between 
different ground states. 



region where J is too small or large so that it would not 
be important to discuss the properties of L1V2O4. There- 
fore, we do not consider these Y E states in more detail in 
the following. 



We represent the eigenstates by a usual notation 2S+1 T 
and if necessary, we will also write the electron number 
explicitly as 2S + 1 r™ d . The five phases correspond to l E, 
7 Ai, 1 A\ : 3 7i and another 1 E. On increasing \t^\, 1 A\ 
state disappears. The other states seem to be robust 
against the variation of t n . 

A wide range of the phase space is covered by 7 A\ 
and 3 7i states. The 7 Ai state has a fully-polarized spin 
moment S = 3 which arises from ferromagnetic double- 
exchange interactions. The 3 Ti state is stabilized in 
the competition between double-exchange and antiferro- 
magnetic super-exchange interactions. Electron configu- 
rations of these ground states are schematically depicted 
in Fig. [6] As is easily seen in (c), the 7 Ai phase cor- 
responds to the case of "strong" Hund's coupling; J is 
larger than the level separations. As shown in (d), the 
3 Ti phase corresponds to the case of "moderate" Hund's 
coupling, since the Hund's coupling is effective only in 

the ' molecular orbitals. There appear two X E states 
in the phase diagram. These states are located at the 



To compare with the = 6 case, we also show the 
ground state phase diagram for rid — 5 and 7, in Fig. [7] 
As is seen, ^ d=5 and 2 T 2 " d=5 states, and 6 E nd = 7 and 
4 A™ d = 7 states are the ground states in a wide range of the 
parameters. It is illuminating to notice that the A A^ d ~ 7 
state is obtained by adding an electron to the molec- 
ular orbital in the 3 T" d=6 state and the 2 T% d=b state is 
obtained by removing an electron from the T 2 ( } orbital. 
This is easily understood in Fig. [BJ In a similar way, 
6£n d =7 s t a t e i s obtained by adding an electron to the E 
molecular orbital in the 7 A" d=6 state and 6 A™ d=5 state 
is obtained by removing an electron from the Ai orbital. 
The fully polarized state a T™ d ~ 7 does not appear in Fig. 
[7] (c). With increasing |t^|, the energy of 8 T" d=7 state 
decreases around the area U > 3.0 eV and J ~ 1.7 eV 
and then becomes the ground state as shown in Fig. [7] 
(d). 2 T™ d-7 state corresponds to the state obtained by 
adding a molecular orbital electron to 1 E rid=G . As 
in the cases of 1 E ni=6 , we will not further discuss the 
region where J is large. 
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FIG. 8: Low-energy level scheme and the quantum numbers 
2S+1 r in each n d . U = 1.5 eV, J = 0.2 eV and U = -0.085 
eV. Arrows shown by solid (dotted) line represent T^~^ (Ai) 
one-particle excitations. For simplicity, we show only the ex- 
citations from the ground states for each n d . 



C. Low-energy spectrum of one tetrahedron 

In this subsection, we proceed to detailed investigation 
of the electron configurations of low-energy states. We 
will justify the schematic picture of Fig. [5] by checking 
whether the same picture applies to the ground states in 
other rid subspaces. We will concentrate on the 3 T 1 IId_6 
ground states since as will be shown in Sec. IIIID1 
the high-temperature thermodynamic properties of this 
phase are similar to the experimental data of L1V2O4. 

The low-energy part of Figs. 0] (a)-(c) is schematically 
shown in Fig. [51 Not shown in Fig. |H the ground state 
of nd = 4 is 1 A\ and the energy gap to the first excited 
states is large, 0.253 eV. Considering this and quantum 
numbers, the ground state of rid = 4 can be considered 
as a "closed shell" state which corresponds to the fully 
occupied E orbitals in the sense of the schematic picture 
in Fig. [5] (a). Starting from this, the low-energy spec- 
tra for rid > 5 can be successively constructed by adding 
electrons in the T 2 ' or A\ molecular orbital as depicted 
in Fig. [5] Note that all the ground states are constructed 
by adding electrons in the ' orbital. Indeed, the fol- 
lowing group theoretical arguments justify this picture. 

In order to characterize the low-energy spectrum in 
more detail, it is important to identify the quantum num- 
bers of the ground states in each rid subspace shown in 
Fig. [SJ As for the total spin, it increases by 1/2 upon 
adding one electron. This simply means that electrons in 
the orbitals tend to align their spins, i.e., the Hund's 
rule. The symmetry of the orbital part can be also un- 
derstood by starting from ground state in the rid = 4 
subspace. This has the closed-shell electron configura- 
tion and therefore the symmetry of Ai-representation. 
The ground state in the rid — 5 subspace is constructed 
by adding one electron in the ^ orbitals. Its symme- 
try is given by the product of two representations, one 
for the starting many-body wavefunction and the other 
for the molecular orbital of added electron. This case is 



simple and the result is A\ ® T 2 = T 2 . 

Next, the ground state for rid = 6 is constructed simi- 
larly by adding the second electron in the T^ - ' orbitals to 
the ground state T!^ -5 . Thus, considering the decompo- 
sition of product representation T 2 <S>T 2 = Ai®E®Ti®T 2 , 
we have four possibilities of orbital symmetry for the 
ground state in rid = 6 space. However, the spin part 
is triplet (S = 1) due to the Hund's rule and its wave- 
function is symmetric, and therefore we should choose an 
antisymmetric representation in T 2 ®T 2 . This is indeed 
unique and T\. Thus, the wavefunction of rid = 6 with 
5=1 should have a T\ symmetry if the ^ orbital plays 
a role of the one-particle excitations. The other states, 
X A\, l E and 1 T 2 appear as excited states as shown in Fig. 

Finally, the ground state for rid — 7 is, once again, 
constructed by adding the third electron in the T^"' 
orbitals. Since the ground-state wavefunction has spin 
3/2 as predicted by the Hund's rule, all the three 
orbitals are occupied by electrons, aligning their spins. 
The symmetry of this state is 4 A 2 . In the group theo- 
retical language, this A 2 is understood from the relation 
Tx T 2 = A 2 © E ®Ti @T 2 . The states other than A 2 
also appear as excited states as shown in Fig. [5J 



D. Spin susceptibility and entropy 

Experiments of magnetic susceptibility indicate that, 
at around T ~ 500 K, the effective moment of vanadium 
ion S 2 S changes from ~ 1.5-1.75 at high temperatures 
to ^0.9 at low temperatures^. This behavior was inter- 
preted as Kondo like screening by Hopkinson et al.^ but 
it is important to check whether alternative explanations 
are possible. The Weiss temperature 8 also changes at 

~ 500 K 8 ' 9 ; it is estimated as 8 500 K by fitting the 

results at 600-1000 K and 8' 30 K at 80-300 K. It is 

valuable to calculate the temperature dependence of sus- 
ceptibility from our one-tetrahedron data and compare 
the results to the experimental data. 

The spin susceptibility Xs(T) per vanadium site is 
given by X s(T) = g 2 li 2 b (S 2 ) /3N S T with g, n B , N„ the 
electron's g-factor, Bohr magneton and the number of 
sites, respectively. Here, (• • ■) denotes grand canonical 
average at temperature T with keeping the average elec- 
tron number at 1.5 per site, and we evaluate this by av- 
eraging over the truncated Hilbert space of rid = 5, 6 and 
7 as explained in Sec. MI Al In order to obtain Xs(T), 
we calculate {S 2 } with varying temperature. 

Figure [9] shows the temperature dependence of the 
inverse spin susceptibility x~ (T) for several J values. 
In the high-temperature regime, the Xs(T) follows a 
Curie- Weiss law irrespective of J values. There is a 
clear crossover marked by arrow at the temperature 
Icross, where the slope of xirH^ 1 ) changes. The high- 
temperature Weiss temperatures 8 is estimated in the 
region of T > 800 K and the low-temperature Weiss tern- 
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FIG. 9: (Color online) Temperature dependence of inverse 
spin susceptibility. U = 1.5 eV, t w — —0.085 eV and U' = 
U — J. The arrows indicate cross-over temperatures T C ross. 



FIG. 10: (Color online) Temperature dependence of charge 
susceptibility. U = 1.5 eV, t n = -0.085 eV and U' = U - J. 
Inset: J dependence of R = T max /r cr0 ss. 



perature 8' is estimated in the region of 50-200 K. The 
results are listed in Table [TT1 together with the magnitude 
of the effective moment per site S 2 S (S'^). (&') and 
(S 1 ^) are estimated by using the following form: 



o2 
°cff 



Xs(T) 



3(T- 0)' 



(2) 



The magnitude of high-temperature S 2 S is comparable 
to (3/4 + 2)/2 = 1.375 which corresponds to that of the 
atomic mixed- valence limit with the Hund's rule, d 1 with 
the spin s = 1/2 and d 2 with s = 1. It is noted that 
the high-temperature S 2 S has only weak J dependence. 
Xs(T) for J = 0.3 eV is qualitatively in good agreement 
with the experimental results^ in three points: 0, T cross 
and S 2 S . The experimental values are O ~ —500 K, 
Tcross ~ 500 K and S 2 S ~ 1.5-1.75^. The deviations 
of S 2 S would be due to the ferromagnetic correlations 
beyond the present calculations. |0'| and S'Jg are smaller 
than high-temperature |0| and S 2 S , respectively. The 
values of S® s reflect the ground and low-lying excited 
states. For J = 0.2, 0.3 and 0.6 eV, S' 2 ^ is very close to 
that for the ground states. Although the ground states 
are spin singlet for J = 0.0 and 0.4 eV, the low-energy 
excited states with spin S ^ contribute S 2 S > in the 
temperature range where we fit the data. The \s (T) for 



TABLE II: Curie- Weiss temperature and effective moments. 
& and iSgfj are estimated by fitting data in Fig. [5] at T > 800 
K. Except J = 0.4 eV, 6' and S'^ are estimated by fitting 
data in Fig. H] at 50-200 K. 9' and S' e % for J = 0.4 eV are 
estimated at 90-160 K. 
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FIG. 11: (Color online) Entropy per tetrahedron unit vs tem- 
perature. U = 1.5 eV, U = -0.085 eV and U' = U-J. 



J = 0.6 eV has a clear ferromagnetic behavior due to the 
double exchange mechanism, see > in Table [Til The 
behaviors at very low temperatures depend on the total 
spin S of the ground state for each parameter set. The 
low-temperature upturns for J = 0.0 and 0.4 eV reflect 
that the ground states are spin-singlet states. 

We find that the cross-over temperature marked by 
arrow in Fig. [5] is related to the energy scale of charge 
fluctuations. The charge susceptibility Xc(T) = (n 2 d — 
(n-d) 2 )/T is shown in Fig. [TO] The peak position of 
Xc(T) (= T max ) is related to T cross of spin susceptibil- 
ity in Fig. [5] The inset of Fig. [TO] shows the ratio 
R = T max /T cross . Apart from large J region, R is nearly 
constant and R ~ 3. From this, we can understand that 
the crossover in spin susceptibility arises from charge fluc- 
tuations at least from small to moderate J values. This 
interpretation is consistent with the effective moments in 
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Table Ulland suggested by the early exact diagonalization 
study 29 • The interpretation of the crossover attributed to 
charge fluctuations in a tetrahedron is valid at least when 
tetrahedron coupling is weak, e.g., at high temperature. 

Figure QT] shows the entropy S(T) per tetrahedron as 
a function of temperature. The finite values at zero- 
temperature are due to the degeneracy of ground states 
in rid = 6 subspace. The experimental data of entropy 
at 100 K from the specific heat data is ~ 5fcslog2 (ks- 
Boltzmann constant) per tetrahedron,A and this is larger 
than the present results at 100 K. Since we ignore inter- 
tetrahedron correlations in the present calculations, it 
is not adequate to discuss the low-temperature entropy 
quantitatively. We note that there still remains large en- 
tropy (more than k B log 9 at 100 K) for J < 0.4 eV. This 
low-energy entropy might become an origin for heavy 
fermion behaviors in L1V2O4. This point will be dis- 
cussed in Sec. HVBI based on an effective model for cou- 
pled tetrahedra. 

IV. LOW-ENERGY EFFECTIVE MODEL OF 
ONE TETRAHEDRON UNIT 

In order to discuss low-energy properties of LiV204, 
one has to notice that among several ground states of 
one tetrahedron unit, 3 Ti phase has both spin and orbital 
degrees of freedom. The magnetic susceptibility and the 
entropy calculated for this phase capture the character of 
the experimental results at high temperature. Therefore, 
we now focus on the 3 Ti phase and discuss its low-energy 
properties in detail. To describe metallic behaviors of 
L1V2O4, it is important to examine one-particle excita- 
tions in this phase. We will construct an effective Hamil- 
tonian for one tetrahedron unit and demonstrate that the 
T 2 '-orbital electrons only are sufficient to describe the 
low-energy one-particle excitations of Hamiltonian ([T]) in 
3 Ti phase. This construction can be regarded as a pro- 
cedure of a real-space renormalization group^ 2 - Based on 
the results obtained in this section, we will proceed to 
the next procedure of the renormalization group in Sec. 

13 

A. One-particle excitations 

Let us now investigate one-particle excitations in the 
3 Ti phase in detail. First, we examine which molecular 
orbitals in a tetrahedron play a dominant role in the one- 
particle excitations upon changing electron number n ( i = 
6^7 and rid = 6 — > 5. To this end, we define matrix 
elements A 7 ^, by 

A TT'= (9n d +l\4"t\9n d ){9n d \dT"t\9n d +l)- (3) 

Here, dL with T = E, T^~\ Ai, T Y and T 2 (+) represents 
the d-electron creation operator in the molecular orbital 
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FIG. 12: (Color online) Occupation numbers of molecular 
orbitals in the ground states in each rid for J = 0.3 eV and 
U = -0.12 eV. (a) U = 1.5 eV. (b) U = 3.0 eV. 

basis with the spin a =|. g rid denotes the ground states 
in the rid subspace. App, can be easily calculated for the 
noninteracting case. We calculated Ap r , and A® r ,, for 
the 3 Ti phase at U = 1.5 eV, J = 0.2 eV and t n = -0.12 
eV. The matrix element {gn d +i\d\- a \g nd ) for T = T% is 
about ten times larger in magnitude than the others. We 
find that the largest eigenvalue A max of App, is suppressed 
to about 80 % of the noninteracting value A^j^. 

(n d = 5) A max ~ 1.6 < A^ a c x = 2, (4) 
(n d = 6) A max ~ 0.80 < A££ = 1. (5) 

It is important that the eigenvector of the largest eigen- 
value has almost all weights in T 2 ' components (~ 
99%). Thus, we obtain the short-range contribution 
of renormalization factor Z ~ 0.8 in the present one 
tetrahedron calculation. The result for U — 3.0 eV, 
,/ = 0.2 eV and t„ = —0.12 eV is that the matrix el- 
ement |(gn d +i|drerl.9«d)l f° r ^ = is still about ten 
times larger than others and we obtain Z ~ 0.66. These 
results indicate that the T% orbitals play a dominant 
and important role for one-particle excitations between 
the ground states of each rid subspace. 

We also calculate the electron occupation number of 
each orbital J2 a &r Y,A9n d \d\ a d aa \g n<l ) for each n d sub- 
space, and the result is plotted in Fig. [T21 The value 
is averaged over degenerate ground states. The occu- 
pation number of the T 2 ^ orbital increases by nearly 
one when rid increases by one. This is consistent with 
the analysis of ASp,. As discussed in Sec. MI B[ a sim- 
ple picture of the 3 T" d_6 ground state is the fully oc- 
cupied E orbitals plus partially filled orbitals. The 
result of orbital occupation confirms this picture but also 
shows that a non-negligible number of electrons occupy 
the high-energy one-particle molecular orbitals such as 
Ti and T 2 (+) , suppressing the i?-orbital occupation from 
four. The larger correlations U, the larger number of 
electrons occupy the high-energy molecular orbitals as 
shown in Fig. [T^] (b). Thus, the one-particle excitations 
should be regarded as dressed quasiparticles in the Lan- 
dau's Fermi liquid picture rather than "free" electrons. 
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to (eV) co (eV) 

FIG. 13: One-particle spectrum — ImGr(a>+i<5) vs u. U = 1.5 
eV, J = 0.3 eV, U = -0.12 eV and 8 = 0.001 eV. (a) F — E, 
(b) T 2 (_) , (c) Ai, (d) Ti, and (e) T 2 (+) . (f) Detail structure of 
the low-energy region. The quantum number 2S+1 r for the 
corresponding peak is also shown. 

These investigations are directly checked by calculat- 
ing Green's functions Gr(w + i<5) in the molecular-orbital 
basis, where 5 is infinitesimal constant. Gr(ui + i5) is de- 
fined as a Fourier transform of a retarded Greens function 
G r (t) at T = 

G r (t) = -i8(t) |{d r (i), 4(0)}\g nd ) (n d = 6),(6) 

where 8(t) is Heaviside's step function, {■ • ■} denotes an- 
ticommutator and we omit the spin index a. We show 
the one-particle spectral function — ImGr(w + i5)/ir with 
6 = 0.001 eV in Figs. Ql (a)-(e). Note that the scale 
of the vertical axis is different for each figure. There are 
large peaks in low-energy region for T = E, T 2 and A\. 
On the other hand, there are no large peak in low-energy 
region for T = T\ and T 2 but broad incoherent compo- 
nents in the high-energy region. This also agrees with the 
simple picture in Fig. [6] In Fig. [13] (f), we show the low- 
energy part of total spectral weight — ImGtott^ + id)/Tr 
defined as G t ot(w + iS) = J2r Gr(w + iS). As expected 
from Fig. the lowest-energy excitations are those of 
T 2 orbital and the corresponding peaks are very large. 
Another important point is that the peak of the Ai or- 
bital is also large and located at low energy. 



B. Effective Hamiltonian of one tetrahedron unit 

The results obtained in the previous subsection show 
that T 2 ' orbitals only are sufficient to describe the low- 
energy sector of 3 Xi phase. In this subsection, we con- 
struct an effective Hamiltonian of these T 2 ' electrons, 
and determine its interaction parameters. 

The general Hamiltonian with spin rotation symmetry 
of T 2 ' electrons at the tetrahedron n, retaining only the 
two-body interactions, should be written as 

era a 
a>/3 <j<j' 

+ ^ H a iaT a Li a «/5ian/3T + C, (7) 



where C is a constant and n naa = a' nacr a nacT with a = a.b 
or c (see, Appendix [A)) . a\ iarj creates a "quasiparticle" 

of T 2 ' orbital at a tetrahedron n that is dressed by the 
interactions and its vacuum corresponds to the 1 A^ d ~ . e 
is the one-particle energy level. The interaction parame- 
ters U, U' and J are the molecular-orbital version of the 
coupling in the £ 2 g Hubbard model |T]) and now the pair 
hopping term T is also generally generated. 

In order to check the validity of Hamiltonian ([7]), we 
compare the numerically calculated eigenenergies (£ num ) 
of four-site case of the original Hamiltonian (fTJ) with one- 
"site" eigenenergies of the effective model (|7|). The low- 
energy eigenvalues are listed in Table HTT1 The number in 
the fifth and sixth columns is (H^g et ) / E num , which mea- 
sures the validity of Hamiltonian ([7]). The results are 
very close to unity and the validity of Hamiltonian |[7J is 
quantitatively proved. The estimated values of the inter- 
action parameters turn out to be smaller than the bare d- 
electron interactions by the factor 1/10 ~ 1/5. This is be- 
cause the orbitals are extended over four sites and there is 
a reduction of energy scale by the one-tetrahedron renor- 
malization factor Z as discussed in Sec. IIV AI The pair 
hopping term T is also induced in this effective model 
but its strength is weaker than the others. 

We note that a few states in Table Mil cannot be de- 
f— 1 — 
scribed by only T 2 orbitals, and these states correspond 

to Ai-orbital excitations as shown in Fig. [8] Although 
we can also construct an effective Hamiltonian including 
these A\ orbitals, we do not try to do this, since the 
model will become too complicated. Indeed, this simpli- 
fication is not so bad, since none of the "ignored" states 
in Table HTT1 is the ground state in any space. It should 
be noted that the above argument does not hold near the 
phase boundary. 
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TABLE III: Comparison of the model Hamiltonian ([TJ and 
the results of truncated exact diagonalization for U — 1.5 eV 
and U = -0.085 eV. (H^ et )/E nurn are shown in the fifth and 
the sixth columns for J = 0.2 and 0.3 eV. The five parameters 
indicated and a trivial constant term are estimated by using 
six "input" states. The states with (rid, S, V) = (6, 1,T2) and 
(7, |,Ti) cannot be described by H e g alone. 



rid 


o 
o 


1 


rrltct s~ t 

H cB - C 


J — U.2 ev 


0.3 ev 


4 





Ax 





input 


input 
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1 
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T 2 
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input 


input 


6 


1 


Ti 


2e + U' — J 


input 


input 


6 





Al 


2e + U + 2T 


input 


input 


r> 




1 

1 


±2 




* 
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n 
u 


IT 


Zt -\- u — 1 


input 


input 


6 





T 2 


2e + U' + J 


input 


input 
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Ai 


3e + 3U' -3J 


0.9989 


1.001 
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Ty 


* 
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1 
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3e + 2(7' + U - J + f 


0.9972 


0.9999 


7 


1 

2 


E 


3e + 3(7' 


0.9967 


0.9963 


7 


1 


Ti 


3e + 2U' + U — J — T 


0.9974 


0.9980 


C 


(eV) 






-5.490 


-4.861 


e 


(eV) 






-0.4090 


-0.3332 


t/ 


(eV) 






0.3140 


0.2626 


U' 


(eV) 






0.2989 


0.2502 


J 


(eV) 






0.0215 


0.0257 


f 


(eV) 






-0.0036 


-0.017 



V. EFFECTIVE MODEL OF FOUR 
TETRAHEDRA AND CORRELATIONS OF SPIN 
AND ORBITAL 

In this section, we will construct an effective Hamilto- 
nian describing interacting tetrahedron units in the 3 Ti 
phase. We will then calculate its low-energy eigenstates 
for the unit of four tetrahedra and the spin and orbital 
correlation functions for the ground states. 



A. Effective Hamiltonian for coupled tetrahedra 

In Sec. IIVB1 we have constructed an effective model 
for an isolated tetrahedron unit. We now derive an ef- 
fective model for coupled tetrahedra in the 3 Ti phase 
by including inter-tetrahedron processes. It is inter- 
tetrahedron d-clectron hoppings that couple otherwise 
isolated tetrahedron units. 

In Sec. IIII Al we obtained low-energy eigenstates {|A)} 
in a tetrahedron unit. When tetrahedron units are decou- 
pled, eigenstates of the whole system are simply direct 
products of the tetrahedron eigenstates: | Ai A 2 • • • \n), 
where N is the number of tetrahedron units. The next 
step of the real-space renormalization group procedure 
is to obtain effective couplings between these low-energy 
states. These tetrahedra are coupled by d-electron hop- 



pings between nearest-neighbor pairs of original sites 

tfho P p= E E(*f4 ff ^ + h-c), (8) 

where ((i, j)) indicates that i and j are the nearest- 
neighbor vanadium sites and belong to different unit cells 
(tetrahedra). There, we need the matrix element of elec- 
tron hopping processes in the tetrahedron basis t\ n \ m '- 

= (AriA m |-ffhopp|A^A^„). (9) 

Here, |A„A m ) is a direct product state of two tetrahedra 
n and m. In practice, we need to calculate the matrix 
element of the d-electron creation (annihilation) operator 
AaoidjfSa)- A typical term in Eq. @ is 

(A„ A m \t?f dl aa d il3cr | A^ X' m ) 
= tf(X m \(X n \dl a d ]l3<7 \X' n )\X' m ) 
= tfi-lf^HXnldljX^iX^d^lX'J. (10) 

where i (j) belongs to the tetrahedron n(m) and Py^ 
is the electron number in |A^). Since the matrix ele- 
ments (A„|4 aCT |A^) and (X m \d^ a \X' m ) in Eq. dTUJ) are 
local quantities, we can evaluate them for the wavefunc- 
tions obtained in Sec. IIIII Using t^"^ m obtained in this 
way, we can write our effective Hamiltonian H e f{ as 

Heff = E eA l^")(^»l 
ii x 

+ E E &\X n X m )(X' n X'J. (11) 

(n,m) Xn\mX' n X' m 

Here, 53 (nm) ^ s ^ ne summation over nearest- neighbor 
pairs of tetrahedra and is the energy eigenvalue for 
one tetrahedron which is independent on n. 

In the actual calculations, we take not only the T% ' 
orbitals related to the one tetrahedron effective Hamil- 
tonian ([7]) but also other orbitals such as A\. This gives 
corrections to Eq. ([7]). Later in Sec. IVII1 we will fur- 
ther simplify this effective model (JTTJ) to a more physical 
form. Since the matrix element {X n \d\ a<J \X' n ) is typically 
of the order of ~ 0.3, and the largest hopping term is 
\t a \ = 0.527 eV in our calculations, the order of magni- 

a' a' — 
tude of |t A " A ™| is estimated as t c g ~ (0.3) 2 x 0.5 = 0.045 

eV, the order of (l/10)|t CT |. This value is relatively 
smaller than the charge excitation energy of one tetra- 
hedron A c ~ 0.1 eV shown in Fig. [TDJ Correspondingly, 
the exchange interaction among tetrahedron units are of 
the order of ^ ff /(2A c ) - (0.05) 2 /(2 x 0.1) = 0.0125 eV. 
This is a new energy scale of the low-energy properties 
of this system. The exchange interactions among tetra- 
hedron units will be discussed in Sec. IVII Indeed, the 
values of the various exchange interactions turn out to 
be less than 0.01 eV. 
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Before starting the detailed analysis of this model, let 
us briefly estimate the number of basis states we need 
to keep for this effective Hamiltonian from the view- 
point of entropy. We are primarily interested in the 
low-temperature behaviors of L1V2O4 below the coher- 
ence temperature T* ~ 30 K. For example the entropy 
at around 100 K is S(T ~ 100K) ~ 5k B log 2 ~ k B log 32 
per four vanadium sites determined from the specific heat 
data4 The effective Hamiltonian should have enough de- 
grees of freedom for reproducing this value. 

As discussed in Sec. IIVB1 the average electron density 
implies that the charge subspaces of = 5,6 and 7 are 
dominant local configurations, and it is natural to con- 
sider a few lowest-energy states in each subspace. Here 
we consider the states with partially filled T% ' orbitals 
and count the total entropy per tetrahedron. For these 
configurations, there are 7715 = 6 states in rid — 5 space, 
77i 6 = 15 states in rid — 6 space and 7717 = 20 states in 
rid = 7 space. This restricted Hilbert space corresponds 
to that of Eq. (JT)) in the case of one tetrahedron. Since 
these three subspaces have different numbers of states, 
the probability of finding each charge configuration, P n , 
is not the same to each other but a function of m's: 



P 5 = P 7 = {2 + a)^ =p, 
P 6 = l-2p= {2a- 1 + l)-\ 



(12) 



where a = y m^m-j j m^. The mixing entropy is therefore 
given by 

7 

S m /k B = -^P„logP„ 

n— 5 

= -(l-2p)log(l-2p)-2plogp. (13) 

Adding the contributions from the degeneracy in each 
charge subspace, the total entropy is obtained as 

7 

Stot = S m + k B ^P n log m n 

n— 5 

logm 6 + log(l + 2a) 
= k B log(rn 6 + 2y/m 5 m 7 ) ~ k B log 36.9. (14) 

This value is close to the experimental estimate at around 
100 K4. 

We repeat the same calculation with retaining only 
the ground states in each charge subspace: = 6, 
?7ig = 9 and = 4. This is a minimal set for describ- 
ing charge fluctuations and electron itineracy. Using the 
same formula, we obtain this time S ~ k B log(9 + 4\/6) — 
k B log 18.8. This value is now large enough to repro- 
duce the value at the coherence temperature T = T*: 
~ 2.5fc_Blog2 ~ fcslog5.66. This suggests constructing 
minimal low-energy effective model defined in this re- 
stricted Hilbert space. We can expect that this describes 
low-energy heavy fermion behaviors. We will propose 
such a t-J like effective model later in Sec. IVIII 



B. Ground state of four tetrahedron units 

Now we investigate the ground state when four tetra- 
hedra in 3 Ti phase are coupled by electron hoppings. We 
are primarily interested in the case of 24 electrons in 
the four tetrahedra in total. To this end, we employ 
an exact diagonalization method for the effective model 
(jlip . Because of memory limit of our computer, we cut 
off the high energy states in our diagonalization. We 
first retain states with rid = 5,6, and 7 in each tetrahe- 
dron, which are minimal states to describe charge fluc- 
tuations and thus the inter-tetrahedron superexchange 
interactions. Secondly, in each rid subspace, we retain 
several lowest-energy states only. The number of the re- 
tained states N™* t is tuned depending on the parameters 
in the model, and typically iVjut — 8, N^ ut = 24 and 
iV(T ut = 32. Inter-tetrahedron correlations are included in 
the energy level e\ and wavefunctions {|A)}. The ground- 
state wavefunction is to be obtained as 

\9) = E'E'E'E'n a 2 a 3 a 4 I Ai A 2 A 3 A 4 ) (15) 

Al ^2 A3 A4 

where denotes the sum over the space restricted by 

The approximation of truncating high-energy states is 
controlled by varying N"^ t and we have checked our re- 
sults by increasing the number of retained states. 

First, we show the ground state phase diagram for 
the case of 24 electrons in four tetrahedra. We calcu- 



^c 6 ut 



24 and 



late ground states by setting A^ ut = 8, 
-/VJut = 32 for J < 0.3 eV and the determined phase 
diagram is shown in Fig. Q3] There are three phases: 
l E, 3 Ti and 7 A2- Once again, a state with total spin S 
belonging to T representation of the Td point group is 
denoted by 2S+1 r. The ground states change from mag- 
netic to non-magnetic one as \t v \ increases. This point 
will be explained in Sec. I VII by estimating the exchange 
interactions between tetrahedron units. We should note 
that the phase boundaries do not converge yet with in- 
creasing the cut-off numbers {iV^}. This phase diagram 
shows approximate, rather than precise, locations of level 
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FIG. 14: Phase diagram for four tetrahedra with rid 
and U = 1.5 eV. iV c 5 ut = 8, 7V c 6 ut = 24 and iV c 7 ut = 32. 
retained states are those listed in Fig. [5] 
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crossing. However, we can learn a few important charac- 
ters of the ground state of the four coupled tetrahedra. 
The first point is that in the shown region of the t^-J 
parameter space, these three states, X E, 3 T\, and 5 A 2 , are 
the three lowest multiplets and their energy separations 
are very small. The second point is that the tendency 
that the states with large spin appear at the small \t n \ 
region is robust among the different truncation numbers 
used. 

Figure Q~5] shows the dependence of the energy 
of the three states appearing in the phase diagram 
relative to that of 3 Ti for different sets of 7V™ u d t , (a) 
(iV c 5 ut ,iV c 6 ut ,7V c 7 ut ) = (8,24,32), (b) (20,34,44) and (c) 
(32, 34, 44). The parameter set (a) is same as that used in 
Fig. M The ground state for (b) is X E for U = -0.12 eV, 
3 Ti for U = -0.11 and -0.10 eV, and 5 A 2 for t„ = -0.09 
eV, and the ground state for (c) is X E for < —0.10 eV 
and 3 Ti for t n = —0.09 eV. The energy of b A 2 strongly 
depends on t„ compared to that of l E and 3 T\. In the 
region of large \t n \, states with a large spin are energet- 
ically unfavored. Superexchange via higher-energy vir- 
tual states are also present, and some of them generate 
antifcrromagnetic correlations. This is understood by ob- 
serving the region of the 5 A 2 state shifts to the small \tir\ 
region as we increase {-/V cu t}. It is noted that the 3 Ti 
and X E states are almost degenerate in a wide range of 
parameters and it is not conclusive which is the ground 
state within the present calculations. The energy differ- 
ence between the ground state and the first excited states 
is typically 10~ 4 -10 -3 eV. This might mean the existence 
of very low- lying excited states in the limit of large N™* t . 



x10 



> 4 

03 



LU 2 



LU 





-2 



(a) (b) (c) 




~i 1 1 1 1 1 1 1 1 1 1 1 1 

-0.12 -0.11 -0.10 -0.09 

t„ (eV) 



FIG. 15: (Color online) t w dependence of the ground state 
energy with T symmetry Er relative to that of 3 Ti for four 
tetrahedra calculations for U = 1.5 eV and J = 0.3 eV. 
The results are obtained by three sets of truncation schemes, 
(a) (A^ ut ,iV c 6 ut ,iV c 7 ut ) = (8,24,32), (b) (20,34,44) and (c) 



(32,34,44). For U = -0.12 eV in (a) (N* ut , iV c 6 ut , iV c 7 u1 
(8, 34, 32) is used because of additional near degeneracy. 



The matrix element (g\S?S?\g) is calculated by inserting 
SaiA2A3A 4 I A i A 2A3A4)(AiA2A3A4| between 5f and 5? as 
a usual procedure. Then, we calculate (g\S? S?\g) from 
one-tetrahedron matrix elements (A^|5f|A n ) for i S n 
and the wavefunction of the ground state \g). Spin cor- 
relation between the two sites in different tetrahedra and 
that in the same tetrahedron are given as follows: 



C. Short range correlations 

Next, we calculate the spin-spin correlation function 
5(q) for the sixteen sites in the unit of four tetrahedra. 
Here, 5(q) is an equal-time correlation, i.e., a frequency 
integrated quantity, and defined by 



(g\s?s!\ 9 ) = E'E'E'n^n 



S(q) 



N 



No 



1 (g\S?S?\g) 
us 

IE 



exp(iq- (xi-xj)), (16) 
(17) 



where \g) and N g means the index and degeneracy of 
ground states, respectively. N s is the number of lattice 
sites (N s = 16 in the present case) and Xi is the position 
of site i. Note that S? is the spin operator not of a 
tetrahedron unit but at the vanadium site i. Since the 
sum over the ground state degeneracy also includes the 
spin multiplet, the correlation is of the scalar part of two 
spin product: 



J2i9\S?S?\g) = l-Y,(9\Si-Sj\g)- 



(18) 



A2 A' \' 



A1A2 A3A4 A' A' 



x (A 3 |S'f|A 3 )(A4|S , /|A;) forie3andje4, 

(19) 

(g\S?S?\g) = E 'J2' W Ux s xWUx 3K 

A1A2A3 A4A^, 

x (A 4 |S'fS , j 2 |A 4 ) for i e 4 and j G 4. (20) 

Here, we have taken the wavefunction W? x x x as real. 

We show 5(q) for U = 1.5 eV and J = 3 o!2 eV in 
Fig. [TUfor the three different phases. This is calculated 
with the cutoff numbers (iV c 5 ut , iV c 6 ut , iV c 7 ut ) = (8,24,32). 
Since 5(0) is proportional to the ground state expecta- 
tion value of Q2i ^i) 2 , There are notable differences near 
q = for different ground states. There is also difference 
in the spatial anisotropy in large q. In the present calcu- 
lations, 5(q) monotonically increases from the zone cen- 
ter to the zone boundary in all the three ground states. 
This behavior is different from the finite |Q*| ~ 0.6 A -1 
spin fluctuation observed in the neutron experiment.- 
This might be due to the fact that the present 5(q) is 
a frequency integrated quantity, while the neutron ex- 
periment observed a low-energy part of spin fluctuations 
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(a) 1 E -0.12 

(b) 3 Ti -0.11 

(c) % -0.09 

J=0.2 eV, qz=0 



1.0 2.0 

qx (1/X) 



FIG. 16: (Color online) Spin-spin correlation function S(q) 
on (q x ,q y ,0) plane. U = 1.5 eV and J = 0.2 eV. (a) l E 
state (U = ~0.12 eV ). (b) 3 7\ state (U = -0.11 eV ). (c) 
5 A 2 {U = -0.09 eV ). The first Brillouin zone is indicated 
by thick lines. The vanadium-vanadium distance is set as 



= 2.85 A. 



(0.2-0.8 meV)^ 

We also calculate the orbital correlations. The orbital- 
orbital correlation function 5"^(q) is defined by 



S?"<q) = 



(9\OUi)Op(j)\g) 

N a 



(21) 



x exp(iq • (x ; - xj)), 

where O a are orbital operators defined by 0^ a = 
i{d\ za d zxa - h.c.)/2, O ib = i(d\ xa d xycr - h.c.)/2, 4c = 
i( d iya d yz<T ~ h.c.)/2, 5a = (d\ za d zxa + h.c.)/2, 5b = 
{d\ xa d xya + h.c.)/2, 5c = (dl ya d yza + h.c.)/2, 3a = 
(2dl ya d xy(7 - d\ Z(J d yza - di xr7 d zxa )/Vl2, and 3 b = 
{d yza d yza — d zxc7 d ZX(7 )/2 (here the vanadium site i and 
a summation are not shown explicitly). The evaluation 
of the matrix element (g\Ol, (i)Op(j)\g) is similar to the 
case of 5(q). We show for X E ground states the real space 
orbital correlations AI aa (ij) = 52 g (g\Ol(i)OpQ)\g)/N a 
with i = 2 fixed and its Fourier transform S" a (q) in 
Figs. Q7| (a) and (c), respectively. Note that M 4b4b (2j) 
and M 5b5b (2j) is identical to M 4c4c (2j') and Af 5c5c (2j'), 
respectively, where j' is the mirror image point of j with 
respect to (110) plane, and therefore we do not plot the 
latter. As we can see in Fig. [T7] (a), inter-tetrahedron 
correlations are strong for the Os , 0$b and 0§ c compo- 
nents. This is clearly seen as a difference in the average 
of |M Qt *(2j)| for 5 < j < 16 as shown in Fig. Q7] (b). As 
for the wavevector dependence, S 5a5a (q) has a peak cor- 
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FIG. 17: (Color online) Orbital-orbital correlation function in 
the 1 ^ state. t n = -0.12 eV, J = 0.2 eV and U = 1.5 eV. (a) 
Real space orbital-orbital correlation functions M at *(2j) be- 
tween site 2 and another site j . The site indices are indicated 
in (d). (b) Average values of inter-tetrahedron correlations 
defined by £)£ 8 |M a< *(2j)|/12. (c) S D 5a5a (q) on (q x ,q y ,0) 
plane. 



respondingly at q = as shown in Fig. [T7](c). The other 
modes of orbital fluctuations have similar q-dependence 
within the first Brillouin zone. We find similar M aa (2j) 
for other ground states and the values of correlations co- 
incide with each other in less than five percents. This 
means that 1 E, 3 Ti and 5 A2 states have very similar or- 
bital fluctuations but the spin correlations are different 
as shown in Fig. 1161 



VI. EXCHANGE INTERACTION BETWEEN 
TETRAHEDRA: SPIN-ORBITAL MODEL 

In this section, we will carry out the second order per- 
turbation calculations in the hopping terms and derive a 
model of Kugel-Khomskii type^ for the spin S — 1 and 
the orbital triplet (T =T±) degrees of freedom, in order to 
investigate the ground states in more detail. The phase 
diagram of four tetrahedra obtained in Sec. IV Bl will be 
explained in terms of various exchange interactions such 
as pure magnetic, pure orbital, and coupled magnetic and 
orbital exchange interactions. Characteristic orbital con- 
figurations coupled to spin degrees of freedom in terms 
of tetrahedron units will be discussed in the final part of 
this section. 
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A. Exchange Hamiltonian 

In order to uncover the obtained ground states and 
their properties, we investigate various exchange inter- 
actions of spin and orbital degrees of freedom between 
different tetrahedra. To examine orbital and spin corre- 
lations, we temporarily neglect charge fluctuations and 
consider 3 Ti multiplet in rid = 6 space at each tetrahe- 
dron. As in a usual manner, we carry out a calculation 
of the second order perturbation in the hopping terms 
(the last term in Hamiltonian (fTTj) ). and derive Kugel- 
Khomskii type exchange interactions 33 of spin and orbital 
degrees of freedom. In the second order perturbations, 
nine states of 3 Ti multiplet are used as initial and final 
states, while eight states in rid = 5 and thirty-two states 
in rid — 7 are kept as virtual states. For rid = 5 and 7 
states, we keep states with 5 = 1/2 and 3/2 in the low- 
energy spectra, since the unperturbed states are those 
with S — 1 for rid = 6. 

We assign the state whose orbital is on the plane^ in- 
cluding the bond (n-m) as T z = one (= |0)). The other 
two states are assigned to T z = ± (= |±)). We show 
in Fig. [18] (a) an example of this assignment. We use 
simplified notations for the orbital label in T\ represen- 
tations hereafter such asxy = (xy + ciz)(x 2 — y 2 ) and so 
on. It should be noted that the same orbital is assigned 
to different Testates depending on bond directions, as 
depicted. We use for the orbital part eight operators, Ij, 
(/j = 1,2, and 8). For orbital degrees of freedom, 
we introduce a representation that depends on the bond 
direction. Let us consider a bond and orbitals at the 
ends of it. /i = 1, 2 and 3 correspond to the pseudospin-1 
operator T Xl T y , and T z , respectively. For /i > 4, we 
define T 4 = {T x , T y }, T 5 = {T y ,T z }, T 6 = {T Z ,T X }, 
T 7 = T 2 - T 2 , and T 8 = [2T 2 - T 2 - T 2 )/V3. For the 
spin part, we use standard spin-1 operators S a (a = x, y, 
and z). Using these operators, the exchange Hamiltonian 
between tetrahedra n and m (bond (n-m)) reads 

HT = E {[Pr(nm) + Jr(nm)S(n)-S(m) 

x T^{n)T u (m)}, (22) 

where T^n) (S(n)) means the orbital (spin) operator at 
tetrahedron n and To (n) = y / 2/3. The pre-factor of Jf 
is just the normalization. Jj 50 is nothing but the origin 
of energy and we set J"° = 0. 

Due to the symmetry of the T\ orbital, selection rules 
exist for J^ v and J^ v and some elements vanish. There 
are two types of symmetry operations which are used to 
reduce the number of independent coupling constants, 
(i) mirror: |+) «-> |— ) for both n and m sites simul- 
taneously, and (ii) C% rotation: n <-> m. First, un- 
der the operation (i), operators T3, T5 and Tq at each 
tetrahedron change their sign while the others do not. 
Thus, the products including one of the former group, 
for example, T^(n)T-j(m)^ cannot appear in the exchange 




site 1 : (xy+c^z)^-^) = xy 



(a) (b) 

FIG. 18: (Color online) (a) An example of orbital config- 
uration. Wavefunctions of Ti are graphically drawn with 
ci = 0.4. The site indices correspond to those in Fig. [1] 
(b) Thick line with arrow indicates the bond favored by jf 8 . 
The direction of the arrow corresponds to that in the graph 
in Fig. EE] (b) Details of the definition T z for bond (1-3). 

interactions, therefore J 37 (nm) = 0. Secondly, under 
the operation (ii), T p (n) (/i =1,2,5, and 6) are trans- 
formed to — Tfj,(m). The others change their site index 
but do not change their sign. From this, the terms includ- 
ing one of them appear in antisymmetric combination, 
T M (n)T„(m) — T w (n)T^(m). The other terms appear in 
symmetric combination (n)T v (m) +T v (n)T ft {m). This 
leads to = —J"^ for the first case, while J^ v = J" 11 
for the second case. Using these properties and the fact 
that the exchange interactions are real 3 ^, it turns out 
that the number of independent couplings is 13 in J^ v 
and 21 in J^ffi. 

We now calculate J^ v and from the numerically 
obtained Hg™, and to this end we use the following Fierz 
identities, 

Tr[S a {n)S b (n)} = 3 x 28 abl (23) 
Tr[T M (n)T„(n)] = 3 x 25^, (24) 

where Tr is taken over in both spin and orbital spaces of 
one tetrahedron. Using Eqs. ([23]) and ([24]) . we obtain 

Om) = ^^'[^(^(m)^], (25) 

JfN = ^^Ty[S(n)-S(m)T>)T*(m)iCT6) 

where Tr' is taken over in both spin and orbital spaces 
for two tetrahedra n and m. In Fig. 119] we show 
and J for the (1-2) bond as a function of t^. Since 
= ±J^ for a = 1 and 2, we plot only one of 
them. It is found that all couplings are smaller than 
10 meV, which is consistent with experimental results 
for the Weiss temperature 9' ~ —40 — 30 K estimated 
below 400 Ki&2. 

The coupling J® is pure spin exchange and decreases 
with decreasing indicating the enhancement of fer- 
romagnetic processes. We can explain this tendency by 
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-0.13 -0.12 0.11 -0.10 -0.09 
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FIG. 19: Spin and orbital exchange couplings Jf(12) and 
J£"(12) as a function of U. U = 1.5 eV, J = 0.2 eV. The three 
smallest J£"(12) and Ji°(12) are not shown for simplicity. 5 ga 
denotes the ground state spin for each t n . 



coupled four-tetrahedron system, 

l<n<m<4 

and calculate a few lowest-energy states. The result is 
that the ground state is X E state for J = 0.2 eV and 
—0.13 <t„< —0.09 eV. By comparing this result to the 
phase diagram of Fig. 1141 it turns out that the pertur- 
bative calculations underestimate the ferromagnetic ex- 
change coupling as is easily understood by observing the 
lack of double-exchange interactions in Hamiltonian (|2T[) . 
The results of the present perturbative analysis is similar 
to the cut-off scheme (c) in Fig. [T5l except = —0.09 eV 
where the ground state is 3 T\. Then, in order to check 
whether we can explain the phase diagram in Fig. [14] by 
the exchange model (|2"T1) , we carry out the same calcula- 
tion by replacing perturbatively calculated J™ by J^s- 
By introducing an effective pure magnetic exchange in- 
teraction J^ff, we can incorporate enhancement of ferro- 
magnetic correlations. We show the three lowest-energy 
eigenvalues obtained in this way in Fig. [50] The three 
lowest states are indeed those appearing in the phase di- 
agram Fig. 1141 This indicates that the present perturba- 
tive calculations capture the essential part of this system. 
As decreases, the ground state changes from X E to 
5 A2- However, the 3 Ti state does not become the ground 
state with varying J®° only, and therefore we would need 
more complete manipulations of the exchange coupling 
constants. 



examining important virtual processes. When adding 
an electron of A\ orbital to the 3 T™ d_6 ground state, 
we obtain basically 4 T" d ~ 7 state which is the lowest ex- 
cited state in rij = 7 subspace (see Fig. [5] and Table 
HIT|) . Since this virtual state 4 T" d=7 has spin 3/2, 4 T" d=7 
state contributes to ferromagnetic exchange interactions 
in the second order perturbations. The point is that the 
energy of this state decreases as decreases. There- 
fore, the ferromagnetic interactions are enhanced. As 
1^1 increases, the energy of 4 T™ d state increases. In 
the large \t^\ region, the antiferromagnetic exchange in- 
teractions generated via the excited 2 T" d - 7 , 2 E n *= 7 and 
2 T" d=7 states dominate. As a result, J$°(nm)S(n) -S(m) 
(pure magnetic exchange interaction) notably becomes 
strong among others and this is antiferromagnetic cou- 
pling. Thus, the inter-tetrahedron exchange interaction 
depends significantly on the excitation energy of A\ or- 
bital. This also explains the tendency observed in Fig. 
1141 i.e., magnetic phases appear in small \t n \ regions. 



B. Spin-orbital model: four coupled tetrahedron 
units 

With the obtained couplings Jf" and J^ v , we numeri- 
cally diagonalize the spin-orbital exchange model for the 



C. Orbital wavefunction 

Now we investigate in detail the orbital part of low- 
energy eigenstates of the exchange model H cx for the 
coupled four-tetrahedron system. The important point 
is that among the pure orbital interactions {Jf}, 
Jf s (nm)Tg(n)Tg(m) term is always the largest and 

2.944 i , , , , , , , 1 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



FIG. 20: (Color online) The lowest three energy eigenvalues 
vs J^g /J$ - trr = —0.10 eV and the other parameters are the 
same as those in Fig. 1191 



16 



nearly independent of t„ 
is defined as 



1 



Since the orbital operator Tg 



T 8 = ^=[|+><+|-2 



+ I-X-I 



(28) 



this term favors the bond configurations in which one 
orbital lies on the plane including the bond (|0)) and the 
other does on the plane perpendicular to that (|±)). The 
system of four coupled tetrahedra has 30 such states, and 
four out of the six bonds have the favored configurations 
in each of them. We can illustrate these 30 states by 
simple graph representations. Typical graphs are shown 
in Fig. [21] Vertices of the square represent tetrahedra. 
For each bond satisfying the condition above, we draw 
an arrow which ends at the vertex (tetrahedron) where 
the orbital state is local |0). In this representation, there 
is at most one arrow going in a vertex but more than one 
arrows can go out from a vertex. There are two distinct 
types of graphs. The graphs in Fig. (3]] (a) are "closed 
path" graphs and contain two orbitals of Ti multiplet. 
The graphs in Fig. [5T] (b) have a shape similar to lasso 
(rope with a noose at end) and contain three orbitals of 
Ti multiplet. 

The orbital part of the ground states for four tetrahe- 
dra can be well described by linear combinations of these 
30 orbital states. When setting = 0, we can show 
that three lowest-energy orbital eigenstates have Ax, E, 
and T\ symmetries. Once again, symmetry classification 
is useful to understand this. The states of type (a) are 
classified as A\ ®E®T\, and those of the (b) type are 
classified as Ax © A 2 © 2E ffi 3Ti © 3T 2 . In Fig. [H 
we graphically show the basis states of each irreducible 
representation for the type (a). Those for the type (b) 
are shown in Appendix [B] The pure orbital terms J^ v 
hybridize type (a) and (b) states. The states of each rep- 
resentation interact in the Hamiltonian only with those 
of the same representation. Therefore, as far as these 
30 orbital states are concerned, the size of the matrix 
to diagonalize is reduced to, 2, 3, and 4 for Ax, E, and 
Ti representation, respectively, and we can diagonalized 
them analytically. The matrix elements for this restricted 
Hilbert space are calculated in Appendix [Bl It is noted 
that another diagonal interaction J 33 < lifts the degen- 
eracy of type (a) and (b), and favors type (a) configura- 
tions. Taking into account the hybridizations between 
type (a) and (b) states, the representations appearing in 
type (a) would have a lower energy. This explains why 
Ai, E, and T\ orbital states are the three lowest-energy 
states. 



D. Spin-orbital coupled wavefunctions 



orbital states are linear combinations of the type (a) and 

(b) states shown in Fig. [5T] For simplicity, we here 

discuss only the type (a) configurations, since the weight 

of type (a) is about two times larger than type (b) in 

the present parameter sets. The six states of type (a) 

are reduced to three irreducible representations A° rb © 
E orb j-orb of Td point grQup ag shown in Fig Each 

of the three has only one set of basis states and therefore 
these states are automatically eigenstates of any orbital 
Hamiltonian with Td symmetry as far as the type (a) 
states are dominant. When the spin-orbital couplings 
J 2" are switched on, these irreducible representations of 
orbital are to be hybridized to constitute eigenstates of 
the spin-orbital system 7T ex . 

First, we start to discuss X E states. The T<j point group 
symmetry of the system implies that the S = sector of 
spin wavefunctions in four tetrahedra is decomposed to 
two irreducible representations A^™ © i<; s P m as shown in 
Appendix IC 41 Since the ground state considered now 
belongs to E representation, this state is a linear combi- 
nation of E mh © Af in , Af h © E s P in and the E represen- 
tation in E mh © E s ? in = A x © A 2 © E. Our calculation 
shows that, among them, E mh © A s P in and E mh © E s v in 
components are much larger than that of A° rh © E spin . 
These dominant two components are entangled with each 
other, i.e., the wavefunction is not approximated by a sin- 
gle product of spin and orbital parts. This means that 
the spin and orbital are strongly coupled with each other. 

To discuss the correlations of orbital and spin further, 
let us calculate for the ground state the probability that, 
upon fixing the orbital configuration to a given one, the 
two spins S(n) and S(m) have the total spin S nm . Note 
that the other two spins also have the same total spin 
Snm, since the E state is spin singlet. The results for 
two representative orbital configurations are plotted in 
Fig. [53] upon gradually switching on the spin-orbital cou- 
plings. Namely, a control parameter S is introduced to 
replace J 2 v — ► SJ^" with < 5 < 1, and the correlations 
are plotted as a function of S. If there are no correla- 
tions either in orbital or spin parts, the probability is 
1/90 ~ 0.011, and the configurations with much larger 



(b) 



(a) 



L 




xy 
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Now let us go back to the phase diagram Fig. Q3] and 
discuss these three types of ground states. In the previous 
subsection, we have discussed the low-energy orbital part 
in detail. Now we proceed to study the spin part together 
with the orbital one. As shown above, the low-energy 



FIG. 21: (Color online) Typical graphical representations of 
orbital states favored by jf 8 . The site indices are indicated by 
the numbers 1 ~ 4 around the first graph in (a). Examples of 
(a) "closed path" graphs, and (b) "lasso" graphs, (c) Actual 
orbital shapes in two representatives. 
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FIG. 22: Six orbital basis states of type (a) in Fig. [21] classi- 
fied by Td point group symmetry. Normalization factor is not 
shown. 



probability are dominant ones. In the orbital part, each 
type (a) configuration has weight 0.13, while 0.04 for type 
(b) at 5 = 0. Overall difference in weights between Fig. 
[23l (a) and (b) is due to this difference in the orbital 
weights. 

It is noted that spin fluctuations are strongly corre- 
lated with orbital configurations. The position of spin- 
singlet (spin-quintet) tetrahedron pair is correlated with 
local orbital configurations as shown in Fig [23] (a) (Fig. 
l23l (b)). For type (a) graphs, spin-singlet correlations 
are strong in the tetrahedron pair for which the orbital 
energy is not favored, i.e., bonds without arrow in the 
figure. This tendency is understood by noting that, next 
to the largest coupling J 88 , the dominant coupling con- 
stants are J^°, Jf 8 and J^ 8 as seen in Fig. QU Their 
contributions to energy are compared for different spin- 
orbital configurations in Table IIVI The largest coupling 
is the pure magnetic exchange J® which is antiferro- 
magnetic. The others are spin-orbital couplings Jf 8 and 
J° 8 which are both positive. Type (a) states have only 
(0, ±), (±, 0) and ±(1, 1) configurations. The sum of the 
three terms give the lowest energy for spin singlet on 
the bonds ±(1, 1), i.e., there exist singlet correlations be- 
tween bonds without arrow. As for the (0,±) orbital 
sector, the energy balance is more delicate but the max- 
imum spin configuration is stabilized. This is because 
although the pure spin coupling favors the singlet one, 
the energy gain from Jf° and Jf 8 is larger for S\% = 2 
configuration. 

For type (b) graphs, spin-quintet correlations are 
strong on two of six bonds, (1-3) and (2-4), as shown 
in Fig. [23] (d). Since the total spin is singlet, this also 
means that spin-singlet correlations are enhanced on the 
other four bonds as depicted in Fig. [53] (e). Existence 
of ferromagnetic correlations can be explained as follows. 
As we discussed in the case of type (a) graphs, the ferro- 
magnetic correlations are enhanced on the bonds (0, ±) 
and (±,0). Since the total spin is singlet, two quintets 
should not be overlapped. Combining these implies fer- 
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FIG. 23: (Color online) Weight of spin wavefunction on each 
bond vs 5. U = 1.5 eV, J = 0.2 eV and t w = -0.12 eV. 
(a) Case of "closed path" orbital configuration. From the 
symmetry, Note that (S14, S24) is equivalent to (S12, S34). (b) 
Case of "lasso" orbital configuration, (c)-(e) Various spin pair 
correlations. 



romagnetic spin correlations on the (1-3) and antiferro- 
magnetic correlations on all the others. 

The other two ground states, namely, 3 7i and 5 A 2 
states can be understood in the same way. The A2 
state has no component of the spin wavefunction with 
T2 symmetry as shown in Appendix IC 2i proved by the 
symmetry argument. We show in Fig. [M]the spin-spin 
correlation for the three types of ground states. Here, 
instead of the usual spin-spin correlation, we decompose 
it into nine parts each of which corresponds to a different 
orbital configuration on the bond considered. Therefore, 
the plotted value includes the probability of each orbital 
configuration. Summing up over the nine parts leads to 
the ordinary spin-spin correlation. We can see that in 
both of 3 Ti and 5 A2 states, the spin singlet correlations 
are strong at (±,±) orbital configurations as in the X E 
case. In the state with larger total spin, of course, spin- 
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TABLE IV: Matrix elements of four dominant exchange interactions for each configuration of two-tetrahedron units: the 



total spin S12 and orbital (T z (l), T z (2)). (±, ±) is the representative for (+, 
representative for (+,0), ( — ,0), (0,+) and (0,—). 
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FIG. 24: (Color online) Spin-spin correlation function pro- 
jected to four distinct cases of orbital configurations (S(l) ■ 
S(2)P ) for 1 E, 2 Ti and 5 A 2 states. P a is the orbital projec- 
tion operator to (T z (l), T z (2)) = (r, r') with r, r' = ±, and 
0. The parameters used are t n = —0.11 eV, U = 1.5 eV and 
J = 0.3 eV. 



spin correlation generally becomes more ferromagnetic 
(shift towards positive). For l E and 3 Ti, the spin-spin 
correlations are nearly absent at (0, ±) orbital configura- 
tions. This comes from the bond average in the definition 
of the spin-spin correlation in Fig. 1231 and ferromagnetic 
contributions (quintets depicted in Fig. [22]) almost can- 
cel with antiferromagnetic ones (singlet). 



VII. DISCUSSIONS 

Let us reconsider the effective renormalizcd Hamilto- 
nian. It is represented in terms of tetrahedron variables 
as Eq. pip . We first analyze it numerically and further 
discussed its kinetic and interaction terms separately, and 
found several important processes. Now let us assemble 
these pieces to build a full effective Hamiltonian. Since 
the low-energy physics is concerned, electron hoppings 
are constrained, and among them we keep minimal hop- 
ping processes. Namely, the electron number in each 
tetrahedron is limited to nd=5, 6, or 7, and we con- 
sider only the hopping processes between ground states in 
these subspaces. The processes including excited states 
are neglected. It is possible to represent this constraint 
if we consider only T 2 ''-molecular orbital, as discussed 



in Sec. IIII CI The E orbital is fully occupied in the cases 
considered, and we represent this configuration as vac- 
uum. We have shown in Sec. IIIII that the ground states 
for rid— : 5, 6, and 7 have spin S — 1/2, 1 and 3/2, respec- 
tively. This means that constrained electron hoppings 
generate ferromagnetic double exchange interactions. We 
note that, as we discussed in Sec. IV Al the total entropy 
of this restricted Hilbert space is S ~ ks log 18.8 per 
tetrahedron, which is larger than the experimental value 
at the coherence temperature T* , ks log 5. 66. 4 Thus, we 
expect that this restricted Hilbert space has large enough 
degrees of freedom to describe the heavy fermion behav- 
iors at low temperatures. 

The other terms of the effective model are interactions. 
We can use the same spin-orbital exchange Hamiltonian 
(|2"2")) for this part, but need to modify its coupling con- 
stants. The values of couplings J^ v and J 2" shown in 
Fig ll9l were determined by including all the possible hop- 
ping processes in the second-order perturbation. How- 
ever, now that we treat the hopping processes connect- 
ing ground states as real processes in the kinetic term, we 
need to subtract their contributions from the exchange 
coupling constants. 

These arguments show that a key issue is the compe- 
tition between ferro- and antiferro-magnetic interactions 
coupled with orbital degrees of freedom. This was dis- 
cussed in Sec. IVID1 and we come to a conclusion that the 
low-energy effective model is similar to the t-J model of 
high-temperature superconductor; it is defined in terms 
of a localized spin one with orbital triplet, and mobile 
quasiparticle with Tj ' symmetry on the effective lattice 
(f.c.c). The localized spin and orbital degrees of freedom 
are coupled via exchange interactions between nearest 
neighbors. The hopping processes of the mobile quasipar- 
ticles change local spin and orbital configurations. The 
leading terms of the effective Hamiltonian read 



fleff = ^(-McffiV(n) + [/ off iV 2 (n) 

n 

E [E(*»- Pa »-«m^P + ri.c. 



(n,m) q/3ct 

2 - 



J2 ( o^i"> m ) + J^(nm)S{n) ■ S(r 



1 1 V 



(29) 
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where P is the projection operator to the restricted 
Hilbert space, namely, the ground states of rid = 5, 6, 
and 7 spaces. a naa represents the mobile quasiparti- 
cle with the t\ ■* orbital and the spin a at the tetra- 
hedron n. N(n) is the number operator defined as 
N i n ) = J2aa a tiaa a naa- S and T M are the localized spin 
one and orbital triplet operators of n& = 6 space, respec- 
tively, /^off = -e + (U' - J)/2 and U cS = (W - J)/2 are 
the effective chemical potential and Coulomb interaction. 
The hopping of quasiparticle (a naa ) is renormalized to a 
smaller value about ~ 400 K at most by two factors. 
One is the overlap of T% ^ molecular orbital with a t2 g 
atomic orbital on one site, while the other is the renor- 
malization factor of quasiparticle [Z ~ 0.8 for U = 1.5 
eV). Precisely speaking, exchange processes are present 
not only for pairs of rid = 6 configurations but also other 
configurations with different rid, but we consider in the 
model (f2"9"| only the former ones, since they are domi- 
nant. As discussed before the exchange couplings J£ v 
(a = 1 and 2) are slightly different from in Eq. ([2"3]). 
but their effects are essentially the same as before. The 
difference is that virtual processes via the ground states 
of rid — 5 or 7 configurations are now not counted for 
J£ v . For example, ferromagnetic contributions in spin 
exchange couplings are reduced leading to J® > J® . 

As we noted above, there exist competing interactions 
some of which favor magnetic ground states, while the 
others stabilize nonmagnetic states. Moreover, the mag- 
netic interactions are strongly correlated with the orbital 
ones. There are thirty-fold degeneracies in the orbital 
configurations in the case of four coupled tetrahedra. 
Due to spin-orbital couplings these degeneracies are lifted 
and we investigated which pair of orbitals favors ferro- or 
antiferro-magnetic correlations. These competitions are 
controlled particularly by the energy level of Ai molecu- 
lar orbital. This is because ferromagnetic spin exchange 
is generated by virtual hopping processes including A\ 
orbital and its coupling constant is enhanced when A\ 
energy level becomes lower. The four-tetrahedron calcu- 
lations in Sec. IV Bl showed that tetrahedron degrees of 
freedom (spin 1 and orbital T\) are partially screened by 
the exchange interactions, which leads to nonmagnetic l E 
ground states. It is quite likely that the heavy fermion 
behaviors of L1V2O4 stem from these competitions. Low- 
temperature metallic behaviors in L1V2O4 are dominated 
by correlated one-particle excitations. We expect that 
these competing fluctuations in spin and orbital also, in- 
fluence the coherence of electron dynamics and strongly 
renormalize their quasiparticle weight. A part of the 
renormalization already comes from fast dynamics in the 
tetrahedron unit discussed in Sec. |IV](Z ~ 0.80 (0.66) for 
U = 1.5 (3.0) eV as a tetrahedron unit). It is expected 
that the quasiparticle weight Z is further renormalized 
to a much smaller value when the effects of low-energy 
excitations in the effective model (|29|) are fully taken into 
account. We expect that due to the competing interac- 
tions in (|29p . the low-temperature quasiparticles (if ob- 
tained) are dressed by the spin, orbital and spin-orbital 



interactions and thus become heavy fermions. 



VIII. SUMMARY 

In the following, we review this paper as a summary. In 
this paper, we have investigated the three-orbital Hub- 
bard model on the pyrochlorc lattice in order to study 
the heavy fermion behaviors of L1V2O4. To study which 
type of degrees of freedom plays an important role in 
low-energy dynamics of this model, we have employed an 
approach of real-space renormalization group type. In 
the first stage of coarse graining, block variables are de- 
fined as follows for each primitive unit cell of pyrochlorc 
lattice, i.e., a tetrahedron composed of four vanadium 
atoms. 

First we numerically diagonalized the three-orbital 
Hubbard model and calculated the ground state and low- 
energy excited states in this unit for the cases of electron 
numbers from rid— 4 to 7. The case of rid=Q corresponds 
to the average density in L1V2O4 (d 1 ' 5 per vanadium 
atom), and other cases describe charge excitations. One 
important result is that these low-energy states can be 
represented very precisely by a simple picture of molec- 
ular orbitals. The ground state of the 71^=4 case has a 
closed shell electron configuration of the lowest molecu- 
lar orbital E. The ground states of the rid— 5, 6, and 7 
cases are described as the fully occupied i5-orbitals plus 
partially occupied ''-orbitals in which electron spins 
are polarized due to ferromagnetic Hund coupling. 

Secondly, we derived an effective Hamiltonian for cou- 
pled tetrahedra as for the next stage of the renormaliza- 
tion group procedure. We have performed this, particu- 
larly for the case of 24 electrons in four coupled tetrahe- 
dra, which corresponds to 16 vanadium atoms constitut- 
ing the cubic unit cell of the original pyrochlore lattice. 
This is also a natural choice of unit for block transfor- 
mation in the second stage of the renormalization group 
approach, and we have calculated the ground state and a 
few lowest excited states of the effective Hamiltonian by 
numerical diagonalization. One important result is that 
there appear three types of ground states in a realistic 
region of parameters in the Hamiltonian and also that 
each of them is degenerate either in the orbital sector 
in the spin sector ( 7 Ai) or in both sectors (^Zi). 
It is also important that these three types of states are 
nearly degenerate to each other, and those that are not 
the ground state are the lowest and the second lowest 
excited-state multiplets. 

Thirdly, we examined in detail which processes are im- 
portant for stabilizing these low-energy states in the four 
tetrahedra. There are two types of processes: one is a 
kinetic term and the other is interaction. The former is 
the process of electron hoppings from one tetrahedron 
to another. The interaction processes do not change the 
electron number in each tetrahedron but do change spin 
and/or orbital configurations. We determined the am- 
plitudes of effective electron hopping between a nearest 
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neighbor pair of tetrahedra and found that they are the 
renormalized to a small value, ~ 0.045 eV. Since the ef- 
fective hopping is small, the interactions are short ranged 
in space and the dominant ones are exchange processes 
of spin and orbital degrees of freedom between nearest 
neighbor tetrahedron pairs. In this effective exchange 
process, each tetrahedron is assumed to have six electrons 
and its electron configuration takes one of the degener- 
ate 3 Ti ground states; i.e., three- fold orbital degrees of 
freedom and spin S — 1 remain. Other tetrahedron con- 
figurations are taken into account only as virtual interme- 
diate states of the exchange processes and they are traced 
out. The interaction Hamiltonian consists of pure spin 
exchanges, pure orbital exchanges and also simultaneous 
exchanges of spin and orbital. We used symmetry argu- 
ments to simplify this interaction Hamiltonian and deter- 
mined its form. Spin space is isotoropic in our starting 
microscopic Hamiltonian and therefore the spin exchange 
is Heisenberg type. Orbital space is not isotropic, but 
there are constraints in the orbital exchanges due to the 
symmetries of the lattice and the orbital wavefunctions 
along with the time reversal symmetry. As a result, the 
pure orbital exchanges are simplified to 13 independent 
coupling constants. Including the pure spin exchange 
and spin-orbital couplings, the effective exchange Hamil- 
tonian has 34 coupling constants in total. They are func- 
tions of the microscopic parameters and we numerically 
determined their values by carrying out the second-order 
perturbation in inter-tetrahedron hopping. 

Fourthly, we calculated the ground state and low- 
energy states of the spin-orbital exchange model, par- 
ticularly for the unit of four tetrahedra. We found that 
two sets of special orbital configurations are stabilized by 
the dominant term of the orbital exchange part. They 
are further coupled to each other by subdominant or- 
bital exchange processes to form three low-energy orbital 
multiplets. These three orbital multiplets are also cou- 
pled with spin wavefunctions and form spin-orbital states 
in low-energy region. There, spin-orbital wavefunctions 
are entangled in orbital and spin spaces. This manifests 
strong coupling of spin and orbital degrees of freedom. 
The ground states obtained in this spin-orbital exchange 
model qualitatively agree with those obtained in Sec. 
IV Bl This means that the overall properties of this sys- 
tem are determined by local spin and orbital degrees of 
freedom. 

Finally, combining these results, we have proposed a 
low-energy effective model for LiV204 in Sec. IVIII The 
effective model proposed contains the competitions of 
double- and super-exchange magnetic interactions cou- 
pled with orbital degrees of freedom. Using this effective 
model, we have discussed the origin of heavy fermion 
behaviors in LiV204. To explain heavy fermion behav- 
iors, it is important to identify the origins of large en- 
tropy at low temperatures. In our effective model, the 
entropy arises mainly from the finite spin (S = 1) and 
orbital (triplet) at each tetrahedron of the effective f.c.c. 
lattice. Usually (typically insulating systems with spin 



or orbital moments), these degrees of freedom undergo 
phase transitions. In our effective model, the spin or or- 
bital moments cannot order due to the competitions of 
interactions. In addition to this, the geometrical frus- 
trations in the effective f.c.c. lattice would also suppress 
phase transitions. This means that, after integrating out 
high energy incoherent excitations in the first renormal- 
ization group step, there are still a lot of low-lying inco- 
herent spin and orbital excitations down to low tempera- 
tures and these excitations prevent quasiparticles formed. 
From these, it is expected that the system evolves Fermi 
surfaces and exhibits heavy fermion behaviors below a 
characteristic temperature, at which well-defined quasi- 
particles appear, that would be suppressed by these in- 
teractions. Interestingly, an insulating phase is found 
at high pressure^. This implies that there are compet- 
ing interactions in LiV204 at ambient pressure. It is an 
open question and interesting to explore the microscopic 
aspect of this transition and the relation between the 
heavy fermion behaviors. It is important to analyze the 
low-energy fluctuations in the effective model l|2T))) to see 
whether a heavy fermi liquid state is realized. Elaborate 
large scale simulations are desired for better understand- 
ing of this model and remain as a future problem. 

As an implication of the present approach, we make a 
comment on the temperature dependence of susceptibil- 
ity. In Ref. 28, an independent tetrahedron description 
was applied to fit the susceptibility data at high tem- 
peratures. We can examine this point by calculating the 
energy change of the ground state when four tetrahedra 
are coupled and it is estimated to be ~ 400 K per tetrahe- 
dron. This scale is not larger than the crossover tempera- 
ture of the susceptibility (T cross ~ 500 K for J = 0.2 and 
0.3 eV) and therefore our arguments based on isolated 
tetrahedron remain qualitatively valid, and the crossover 
is mainly due to the suppression of charge fluctuations. 
Of course, inter-tetrahedron spin correlations also con- 
tribute to the temperature dependence of magnetic mo- 
ments and this is also an important future problem. 

We make another comment on the scenarios of the 
Kondo effect or the Mott transition. In these scenarios, 
localized a\ g orbitals play an important role to explain 
the heavy fermion behaviors of L1V2O4. Our result is 
not consistent to such a situation. In the realistic pa- 
rameter space, our calculations show that the density of 
a\ g electron is far below unity per site in the low-energy 
sector. This feature is not consistent with these scenarios 
where the essential point of physics lies in the half filled 
configuration of a\ g orbital. Experimentally, as observed 
by Jonsson et al., L1V2O4 remains a bad metal at high 
temperature^ Moreover there is no signature of logarith- 
mic increase in the resistivity in the whole temperature 
region. These results do not support the Kondo scenario 
in L1V2O4 either. 

In conclusion, we have investigated an effective Hamil- 
tonian of three-orbital Hubbard model on a pyrochlore 
lattice. We have discussed the inter-tetrahedron cor- 
relations and one particle excitations by carrying out 
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two-stage real space renormalization group calculations: 
a tetrahedron unit and then four coupled tetrahcdra. 
We have concentrated on 3 Ti phase of one tetrahedron 
which has spin-one and orbital-triplet ground states. It 
is found that the one-particle excitations in 3 Ti phase 
are described by only ^ molecular orbital even in the 
strongly correlated regime. We have derived an effective 
exchange model in the form of Kugcl-Khomskii model 
with spin one and orbital triplet. Low-energy orbital cor- 
relations are analyzed together with spin-orbital correla- 
tions. It is found that orbital correlations are strongly 
coupled with spin correlations. Finally, we have pro- 
posed an effective Hamiltonian for LiV2C>4 similar to 
a t-J model, in which there are competing ferro- and 
antiferro-magnetic interactions coupled with orbital con- 
figurations together with mobile electrons. These com- 
peting interactions are expected to generate a new small 
energy scale and becomes an origin of heavy quasiparti- 
cles with cooperating with geometrical frustration of the 
pyrochlore lattice. These results would provide a good 
starting point for the further studies of the renormaliza- 
tion group analysis to understand the exotic properties 
in L1V2O4. 
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APPENDIX A: ONE PARTICLE ORBITAL 



In this Appendix we show the wavefunctions for the 
one-particle molecular orbitals. There are twelve states 
as molecular orbitals for one tetrahedron in our model: 
Ai, E, T\ and 2T 2 . Since there are two kinds of T 2 or- 
bitals, these two states can mix with each other. The 
d-electron annihilation operators in the molecular orbital 
basis dr are given as follows (we omit the site and spin 
indices). 



d E „ 
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where the signs are ({a n }) = (+,—,+,—), ({/?„}) = r„ = (a„,/3 n ,7„) coincides with the direction from 
(+,-,-,+) and ({7„}) = (+,+,-,-). Note that 
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FIG. 25: Orbital basis wavefunctions. (a) closed path graphs, 
(b) lasso graphs. See also Fig. I21I (c) examples of ip' = Rip. 

the site n to the center of the tetrahedron. We label 
three states of T\ and two T 2 representations such that 
(Ti a ,T lb ,T lc ) oc ((xy + c 1 z)(x 2 - y 2 ),(yz + cix)(y 2 - 
z 2 ), (zx + c iy )(z 2 - y 2 )), and {t£\ t£\t^) cx ((xy + 



C2 z), (yz + c^'^a;), (za; + c^'y)) for n = 1 and 2, where 

ci and c 2 are constants. The site indices on the right 
hand side of Eqs. (|Alj) - (|A7j) are those in a unit cell and 
indicated in Fig. [lja). 
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APPENDIX B: ORBITAL WAVEFUNCTIONS 



We solve the eigenvalue problem of the orbital exchange 
Hamiltonian in the subspace of these thirty states of type 
(a) and (b). 

The cluster of the coupled four tetrahedra has also a 
tetrahedral symmetry Td and this is useful to simplify 
the eigenvalue problem. As explained in Sec. IVI C\ the 
six states of type (a) are classified to three irreducible 
representations, A\ © E © T\ and they are given as 

\a,Ai) = ■jg(l+R){<p a x+<p ay + <p a ,), (Bl) 
\a,E a ?-yi) = %(l+R)(<p ttx -tp a y), (B2) 
\a,T la ) = ^=(1 - R)<p ax . (B3) 

Here R is the operator that reverses the arrow direction 
in the closed path part, i.e., Rip = ip', and the other 
basis states of the E- and Ti-representations are obtained 
by applying appropriate symmetry operations to these. 
Similarly, 24 states of type (b) are classified to A\ © A 2 © 
2i?ffi3Ti©3T 2 and the representatives of their basis states 
are 
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= ^7|(1 T (¥>13 - ^24 + <P31 - ^42* 



(2) 
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b T w ) 

' ^(3) ? = 775 [(^12 ~ V21 + <^34 - ^43) 



=Fi?(<^14 - <^23 + <^32 - V4l) ■ (B8) 



The other basis states are also generated by applying 
appropriate symmetry operations. 

The orbital exchange Hamiltonian for the four tetra- 
hedra 



In this Appendix, we study the orbital part of the ex- 
change model (|2T|) and explain in detail the calculation of 
low-energy eigenstates in the unit of four coupled tetra- 
hedra. As discussed in Sec. IVI Ci two sets of states are 
favored by the largest term J 88 of the orbital couplings: 6 
states of type (a) and 24 states of type (b). Half of them 
are shown in Fig. [25] (a) and (b) with the arrow repre- 
sentation explained in Fig. [^U Some of these are defined 
with minus sign as a phase factor for later convenience. 
The other half of the states are defined by reversing the 
direction of arrows in the part of closed path and de- 
noted with prime symbol like ip' as shown in Fig. [25] (c). 



-fforb = 



l<n<m<4 



§ Jf " (mn)T^ (m)T v (n) 



(B9) 



has finite matrix elements only between the basis states in 
the same representation. In the subspace of thirty states 
of type (a) and (b), some pairs of coupling constants 
are not independent and it is convenient to introduce 
the parameters K± = (J\ x + jf>) ± (J 22 + Jf ). The 
Jf 8 term gives a constant energy Eq = — 6Jf 8 in this 
subspace. Aside from this constant, the results are the 
following. 
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a. Ai-representation: (dimension 2) 

(a,A 1 \H olh \a,A 1 ) = ±J 33 , (BIO) 

(b, Ai\H olh \b, Ai) = -K+ - 3/3 Jj 18 + aJ\ 7 - ajf - ajf, (Bll) 

(a, Ai\H olh \b, A{) = IK- - 2 7 J 78 + 2/3J x 18 - aJ\ 7 - ajf 2 , (B12) 

where a = \/8/3,/3 = (2/3) 3 / 2 and 7 = 4/V27. 
6. A 2 -representation: (dimension 1) 

(b, A 2 \H olh \b, A 2 ) = \K + - 2 7 J 78 - f3J\ 8 - aJ\ 7 + aJ 36 + aJf . (B13) 

c. E -representation: (dimension 3) 

(a,E\H orb \a,E) = p 33 , (B14) 

(6, E^\H olh \b, E^) = -X+ + I/3J1 8 - \aJ\ 7 + fajf + ±ajf , (B15) 

(6, £ (2) 1^16, £ (2) ) = ±X + + 7J 78 + i/3^ 18 + \aJ 17 - \aJf - \ajf, (B16) 

(a, £|#orb|k, = \K. - 7J 78 + ^Jf + aJ{ 7 - I Jf + aJf, (B17) 

(a, ^|i/ OTb |6, £ (2) ) = ^{K- + 3 7 J 78 - 3f3Jl 8 ), (B18) 

(6, £«|if orb |&, £< 2 >) = ^(jj 78 - (5J\ 8 - aJ\ 7 - aJ 36 + aJt 2 ). (B19) 



d. Ti-representation: (dimension 4) 

4 
3 



a,Ti\H or b\a,Ti) — jJ 33 , 



(B20) 

i^l |-"orb|") - 1 1 / — — 3^+) (B21) 
( 2 )|77 . If, T^( 2 )\ _ Cl ~, 778 1 « 7I8 1 „, 736^ 



6,T 1 (3) |fl orb |6,Ti 



(3)\ 
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(B22) 
3a Ji 2 0B23) 
(B24) 
(B25) 
(B26) 
(B27) 
-aJ 36 ,(B28) 
(B29) 



o.Tilflorblft.T^) = -a(Jl< + Ji ), 

a, T l \H mh \b,T[ 2) )=-lK_, 
a^Ho^T^) = 2 7 J 78 - 2/3 J{ 8 , 

b, T^\H m . h \b,T^) = l 1 J 78 +(3Jl 8 , 

b,TP\Ho*\b,T{ 3) ) = -&J? + Wi 18 

b,TP\Hor b \b,Ti 3) ) = -hJ 7 i 8 -^l S - 
e. Ti-representation: (dimension 3) 



(b 7 T 2 {1) \H olh \b,T 2 

I IT I h 



^4 



(B30) 



(6,T 2 w |if orb |6,^ 

i(^) I IT I k T^) 



iif ++7 J 1 78 -/3J 1 18 -aJ 1 36 , 
(B31) 

(6,T 2 w |i7- orb |6,T 2 w ) = -\K + - aJ\ 7 + ajf, (B32) 

<&,T 2 (1) |tf orb |&,T 2 (2) > = ± 7 J 78 +/3J! 18 , (B33) 

(6, T 2 (1) |H orb |&, T 2 (3) > = | 7 J 78 - aJ 36 , (B34) 

<&,T 2 (2) |i/ orb |&,T 2 (3) ) = -±7 J 78 - PJI 8 . (B35) 

Thus the Hamiltonian is reduced to small matrices and 
the largest size of matrix is four. It is possible to obtain 
analytic expressions of the eigenenergies, but we do not 
write here very lengthy results. 



TABLE V: List of the spin wavefunctions on a tetrahedron 
constructed by four spin-1 states. 



s 


r 


4 


Ai (singlet) 


3 


T 2 (triplet) 


2 


Ai (singlet), ^(doublet), T 2 (triplet) 


1 


Ti (triplet), T 2 (triplet) 





A! (singlet), _E(doublet) 



APPENDIX C: WAVEFUNCTIONS FOR FOUR 
S=l SPINS ON A TETRAHEDRON 



In this appendix, we show the spin wavefunctions on 
a tetrahedron constructed of four spin S — 1. These 
wavefunctions are classified by the total spin 5* and the 
irreducible representation T of point group and listed 
in Table [VJ The point group Td is isomorphic to the 
symmetric group S4 when permutations of tetrahedron 
vertices are concerned, and therefore Young diagrams can 
alternatively be used for irreducible representations, see 
Fig. This is useful particularly when we see the 

symmetries of the wavefunctions. 

In the following, the spin wavefunctions 
are represented by the linear combination of 
|*,(l)s,(2)«,(3)«,(4)>, where s z (n)(= -1,0 and 1) 
represents the eigenvalue for the z-component of the 
spin at the site (tetrahedron) n in Fig. [1] (b). For 
convenience, we write —1 as 1 and list the highest states 
(S z = S) below. 
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FIG. 26: Young diagrams for symmetric group S4. 



1. Wavefunctions for 5 = 3 and 4 

Although we do not discuss the spin wavefunction with 
S = 3 and 4 in this paper, we list the form of the wave- 
functions for completeness. 



S = 4, S z = 4 



%) = |1111). 



(CI) 



5* — 3 . (S'^ 



l 7 r 2n ) =x 



2/./ 



T 2c > = 



|1110) + 1 1101) - |1011) - |0111) , (C2) 
|1110) - 1 1101) + |1011) - |01 11)1 , (C3) 
- |1110) + |1101) + 1 1011) - |0111)j(.C4) 



2. Wavefunctions for S = 2, S z = 2 

The spin wavefunctions for S — 2 are used for dis- 
cussing the 5 A 2 spin-orbital ground states of four tetra- 



hedra in Sec. IvTDl 
V6 



^1) 



\*T, 



2(i/ 



\ 5 T, 



2V7 



1 



1 1 111) + permutations 

1 0011) + permutations , (C5) 
|1010) - |1001) - |0110) + |0101)1 (C6) 
2|1100) + 2|0011) - |1001) - |1010) 

(C7) 

1 1111) + 1 1111) - |im) - |iin) 

|1100) - |0011)j, (C8) 



2V2 

-|0110) - 10101) 
1 



2b) — 



To r ) = -= 



V6 
1 

~V6 
1 

7e 
1 

1 

V6 
1 

76 



1 1111) - 1 1111) + 1 1111) - lllll) 
|1010) - 10101) 



(C9) 



- |1111) + |1111) + |1111) - |1111) 



1 1001) - |0110) 



(C10) 



3. Wavefunctions for S = 1, S z 



The spin wavefunctions for 5 = 1 are used for dis- 
cussing the 3 Ti spin-orbital ground states of four tetra- 
hedra in Sec. IVIDI The wavefunctions are given as 



l 3 Ti 



16 



|3 T; 



? T 2b ) 



|3 T: 



2c 



1 



2y/2 
1 



2V2 
1 



2V2 
1 



2a) — 



|0111) - |1011) + 1 1101) - |1110) + |1110) - 10111) + 1 1011) - |1101) 
10111) - |1101) + 1 1110) - |1011) + |1011) - 10111) + 1 1101) - |1110) 
10111) - |1110) + |1011) - |1101) + |1101) - 10111) + |1110) - |1011) 

1 



10001) + 10010) - 10100) - 1 1000) 



|oiii) + |ioii) - 111011) - |mo) 



1 



2vT0 

1 



|0111) - |1110) + |1011) - |1101) + |0111) - |1101) + |1011) - |1110) 

1 



10100) - 10001) + 10010) - 1 1000) 



10111) + |1101) - 1 1110) - |1011) 



1 



2vT0 

1 



|0111) - |1011) + |1101) - |1110) + |0111) - |1110) + |1101) - |1011) 

1 



10001) + 10010) + 10100) - 1 1000) 



|0111) + |1110) - |1011) - |1101) 



1 



2V10 



10111) - |1101) + |1110) - |1011) + 10111) - |1011) + |1110) - |1101) 



(Cll) 
(C12) 
(C13) 

(C14) 

(C15) 

(C16) 
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4. Wavefunctions for S = 



The spin wavefunctions for S — are used for discussing the E spin-orbital ground states of four tetrahedra in 
Sec. IVIDI The wavefunctions are given as 



-&x 2 — y 2 l 



3\/5 



1 1111) + permutations 



1 1001) + permutations 



sOOs) + \0ss0) + \ssss) - \ssss) - \0s0s) - \s0s0) 



s=±l 



(C17) 
(C18) 



l 1 ^.^) = - \2\ssss) - \ssss) - \ssss) + 2|ss00) + 2|00ss) - |s00s) - |0ss0) - |0s0s) - |s0s0) . (C19) 



*=±i 



When we rewrite these wavefunctions by using direct 
products of two bond spins, e.g., = 3^(1) + S^(2) 
and = S^(3) + <S^(4), we obtain 

|%) = ||22) + ^|00), (C20) 
\ x E x ^ y .) = |11), (C21) 



r^-r») = ^|22)-||00). (C22) 

Here | Si 2 £34) represents the spin singlet state con- 
structed from the bond state with the total spin S12 and 
S34. 
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